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PREFACE 


The  research  described  herein  was  performed  by  The  University 
of  Tennessee  Space  Institute,  Tullahoma,  Tennessee,  for  the  Federal  Avi¬ 
ation  Agency,  Washington,  D.  C. ,  under  Contract  No.  DOT-FA70WA-2260. 

The  Federal  Aviation  Agency  project  monitor  was  Mr.  J.  K.  Power. 

The  work  reported  herein  was  performed  during  the  time  period  November  21, 
1969  through  November  21,  1971. 


This  research  project  encompassed  four  main  areas.  The  principal 
investigator  was  Dr.  B.  H.  Goethert,  Dean  of  The  University  of  Tennessee 
Space  Institute  and  Professor  of  Aerospace  and  Mechanical  Engineering. 
Mr.  Robert  W.  Kanun,  Executive  Assistant  to  Dr.  Goethert,  was  overall 
administrative  coordinator  for  the  project. 

The  various  areas  investigated  and  the  main  investigators  involved 

were: 


Fundamental  Theories  -  Dr.  Y.  S.  Pan,  Associate  Professor 
of  Aerospace  Engineering  and  Dr.  S.  N.  Chaudhuri,  Associate 
Professor  of  Aerospace  Engineering. 

Sonic  Boom  Reduction  by  Focused  Laser  Beam  Techniques  - 
Mr.  Ronald  Kohl,  Assistant  Professor  of  Physics. 

Sonic  Boom  Reduction  by  Cryogenic  Cooling  of  Air  -  Dr.  Marcel 
K.  Newman,  Professor  of  Mechanical  Engineering. 

Engine  Airframe  Integration  with  Special  Emphasis  on  Noncircular 
Engine  Exhausts  and  Jet  flaps  -  Dr.  Heinz  Gruschka,  Associate 
Professor  of  Physics  and  Mr.  Philip  Kessel,  Ph.  D.  Candidate  in 
Aerospace  Engineering. 
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1.  INTRODUCTION 


Even  though  the  U. S.  recently  discontinued  the  SST  project, 
the  problem  of  sonic  boom  attenuation  is  still  of  decisive  im¬ 
portance  because  of  supersonic  flights  of  present  and  future 
military  aircraft  and  future  generations  of  space-flight  vehicles, 
as  well  as  future  commercial  aircraft .  Flight  routes,  particularly 
of  both  domestic  and  foreign  supersonic  commercial  transports 
will  be  restricted  to  over-water  routes  if  the  boom  attenuation 
efforts  are  not  sufficiently  successful;  whereas  the  entire  air 
space,  both  over  land  and  over  water  can  be  opened  for  super¬ 
sonic  air  transportation  if  effective  and  economical  attenuation  . 
of  the  sonic  boom  can  be  achieved.  Thus  the  entire  future  of 
civilian  supersonic  air  transportation  is  at  stake. 

Numerous  theoretical  and  design  studies  as  well  as  extensive 
experiments  have  been  conducted  to  gain  an  understanding  of  the 
physical  phenomena  involved,  and  to  develop  methods  for  pre¬ 
dicting  and  influencing  the  boom  signature  shape  and  intensity 
through  appropriate  aircraft  design  and  operation.  Mainly  during 
the  late  1960’s,  systematic  efforts  with  the  above  objectives 
in  mind  have  gained  momentum  and  much  progress  has  been  made. 

For  instance,  it  is  possible  today  to  predict  satisfactorily 
well  the  sonic  boom  intensity  and  shape  in  the  near,  mid,  and 
far  field  for  current  supersonic  aircraft  designs,  and  to  estimate 
with  adequate  accuracy  the  influence  of  design  changes!  and  of 
environmental  parameters  within  the  framework  of  today’s  con¬ 
ventional  aircraft  designs. 


1.1 


Because  of  the  pressing  urgency  for  finding  an  immediate 
solution  for  th**  sonic  boom  problem,  most  of  the  recent  investi¬ 
gations  and  analysis  methods  were  directed  towards  the  problems 
associated  with  the  planned  first -generation  of  supersonic  trans¬ 
ports.  As  a  result,  unconventional  supersonic  transport , designs , 
employing  novel  configurations  of  volume  displacement,  lifting 
surfaces,  engines  and  exhaust  jets  (jet  flaps)  and  integration  of 
such  components  from  the  viewpoint  of  optimum  sonic  boom  sup¬ 
pression  have  not  received  much  emphasis.  The  University  of 
Tennessee  Space  Institute  and  the  Federal  Aviation  Administration 
recognized,  however,  that  from  a  long-range  viewpoint  the  gat 
in  our  knowledge  concerning  unconventional  configurations  must 
be  filled  by  appropriate  research.  It  was  also  recognized  that 
such  research  would  not  only  guide  the  designers  and  operators 
in  devising  effective  evolutionary  modifications  of  the  first- 
generation  supersonic  aircraft  for  the  purpose  of  better,  noise 
suppression.  It  would  also,  and  probably  more  importantly, 
provide  the  design  base  for  the  next  generation  of  supersonic 
and  hypersonic  transports  for  coping  with  the  sonic  boom  problem 
more  effectively  than  can  be  expected  for  the  first-generation 
supersonic  aircraft ,  even  after  successive  evolutionary  improve¬ 
ments  . 

The  University  of  Tennessee  Space  Institute  therefore  under¬ 
took  a  multiphase  investigation  for  the  Federal  Aviation  Admini¬ 
stration,  the  objectives  of  which  were  defined  in  the  contract  as 


follows: 


"Fundamental  Theories 


This  task  will  be  to  assess  the  numerous  theories  and 
programs  presently  available  with  an  objective  of  defin¬ 
ing  regions  of  applicability  and  areas  of  limitation. 

Summary  descriptions  of  near,  mid  and  far  field  theories 
will  be  compiled  with  statements  on  the  fundamental  physi¬ 
cal  assumptions  and  simple  verifications  of  limiting  factors. 
The  formulation  oi  improved  potential  solutions  for  the 
far  field  will  be  explored,  which  are  currently  visualized 
as  utilizing  not  only  displacement  singularities  in  a 
potential  flow.  Particular  emphasis  will  be  placed  on 
uefining  the  limitations  of  two-dimensional  theory  in 
predicting  sonic  boom  signatures  of  aircraft.  An  explora¬ 
tion  and  evaluation  of  potential  solutions  in  the  vicinity 
of  caustics  will  also  be  made.  To  supplement  the  theoretical 
approach  into  the  right  direction,  some  orientating  wind 
tunnel  tests  will  be  conducted  in  the  later  phases  of  this 
program  area." 

Dr.  Y.  S.  ^n  and  Dr.  S,  N.  Chaudhuri  each  assisted  by  several 

graduate  students,  investigated  separate  phases  of  this  task. 

Details  of  their  findings  are  reported  in  Sections  2.1  and  2.2 

respectively,  of  this  report. 

"Sonic  Boom  Reduction  by  Focused  Laser-Bean;  Techniques 

A  feasibility  study  supported  by  wind  tunnel  experiments 
will  be  conducted  to  determine  whether  focused-laser  beam 
techniques  can  be  used  to  produce  changes  in  the  flow  which 
may  be  relatable  to  equivalent  airframe  shapes.  Realistic 
analyses  of  weight  penalties  associated  with  focused  laser 
apparatus  will  be  included  to  evaluate  the  feasibility  with 
respect  to  aircraft  application." 

Dr.  Ronald  Kohl  conducted  this  feasibility  study,  which  is  reported 

in  Section  2.3  of  this  report. 

"Sonic  Boom  Reduction  by  Cryogenic  Cooling  of  Air 

This  task  will  be  to  assess  briefly  the  quantitative  improve¬ 
ment  potential  and  the  penalties  associated  with  a  radically 
new  boom  suppression  technique  by  making  use  of  air-liquid¬ 
ation  effects.  Disregarding  the  final  practicability  of  this 
novel  scheme,  such  studies  will  stimulate  thinking  along 
unorthodox  lines  and  thus  will  have  significant  indirect 
benefits  which  can  only  be  discovered  by  deviating  greatly 
from  the  path  of  accepted  conventional  designs." 


1.3 


Dr.  Marcel  K.  Newman  determined  during  the  first  period  of  the 
contract  that  the  additional  power  required  to  maintain  large 
external  surfaces  at  cryogenic  temperatures  in  the  Mach  No.  3 
flow  field  was  excessive.  Thus,  with  the  concurrence  cf  the 
FAA  Project  Monitor,  this  work  wan  discontinued  and  is  not 
reported  in  detail. 


It  is  believed,  however,  that  cryogenic  cooling  of  the  air  would 
be  feasible  in  a  configuration  known  as  "space-plane”.  The 
’space-plane"  concept  received  much  attention,  and  extensive 
perfor*:.ance  and  design  studies  were  conducted  in  the  early  1960 
period.  The  propulsion  system  of  the  space-plane  is  of  the  liquid 
air-type  ^.n  which  air  is  liquified  before  it  is  burned  in  the 
combustor.  By  means  of  the  liquif ication  of  the  air,  the  solid 
displacement  of  a  major  part  of  the  airframe  structure  can  be 
eliminated,  and  the  sonic  boom  correspondingly  reduced. 


gEngine-Airlrame  Integration  with  Special  Emphasis  Cn 
Noncircular  Engine  Exhausts  and  Jet  Flaps  - 


This  ..ask  will  be  to  explore  the  potential  and  feasibility 
of  unconventional  aircraft  configurations,  in  which  the 
complete  integration  of  fuselage-wing-Dropulsion  units  is 
approached.  Special  emphasis  will  be  placed  on  examining 
in  epth  the  potential  of  noncircular  engine  exhausts,  pos- 
sibly  m  combination  with  variable  jet  plumes,  for  b'-Ciu  re¬ 
duction.  Besides  other  advantages  of  jet  flaps,  exhausting 
the  processed  air  at  the  trailing  edge  of  the  wing  either 
over  the  entire  wing  span  or  only  over  part  span,  the  reduc¬ 
tion  of  the  rate  of  area  increase  of  the  stream  tube  air 
passing  through  the  engine  may  be  convenient 2y  integrated 
wi„h  the  air  frame  displacement  and  the  equivalent  displace¬ 
ment  due  to  lift.  A  large  number  of  such  highly  integrated 
engine-airframe  configurations  may  become  feasible  especially 
if  stability  considerations  of  the  aircraft  would  be  elimina¬ 
tes,  ana  inste? i  artificial  stability  and  control  relied 


Dr.  B.  H.  Goethert,  principal  investigator  for  this  project,  assisted  by 
Dr.  Heinz  Gruschka  and  Mr.  Philip  Kessel  (A  UTSI  Doctoral  Candidate) 
studied  this  task.  Their  findings  are  reported  in  Section  2.4. 

-Publications  in  professional  journals  based  on  the  research  work  of  this 
contract  were: 


Pan,  Y.  S.  (1970),  "Application  of  Whitham‘s  theory  to  sonic  boom 
in  the  mid-or  near-field",  A1AA  Journal,  Vol.  8,  No.  11, 

pp  2080-2082 

Pan,  Y.  S,  (1971),  "A  method  for  wind  tunnel  investigations  of  sonic 
boom  based  on  large  models",  AIAA  Paper  No.  71-184 

Chaudhuri,  S.  N.  and  Praharaj,  S.  C.  (1971),  "The  near-field  flow 

pattern  of  an  inclined  slender  body  of  revolution",  AIAA  Paper 
No.  71-626. 


^addition  to  these  publications.  Dr.  Pan  presented  Paper  No.  71-184  at 
the  9th  AIAA  Aerospace  Sciences  Meeting,  January  25-27,  1971,  New  York; 
and  Mr.  Praharaj  presented  Paper  No.  71-626  at  the  4th  Fluid  and  Plasma 
Dynamics  Conference  of  the  American  Institute  of  Aeronautics  and  Astronautics 
June  21-23,  1971,  Palo  Alto,  California. 


Additional  publications  and  presentations  are  also  in  preparation. 


2. 


DISCUSSION 


2.1  EXTENSION  OF  THE  CURRENT  SONIC  BOOH  THEORY  AND  DEVELOPMENT 
OF  A  WIND  TUNNEL  TESTING  METHOD  BASED  ON  LARGE  MODELS 

by  Dr.  Y.  S.  Pan  and  UTS I  students  K.  T.  Wang  and  M.  O.  Varner 

2.1.1  INTRODUCTION 

Because  of  the  development  of  supersonic  transports,  the 
sonic  boom  problem  has  been  receiving  considerable  attention 
in  the  past  decade.  Calculations  of . sonic  boom  pressure  signature 
of  a  supersonic  aircraft  have  been  based  mainly  on  Whitham's 
supersonic  projectile  theory  (Whit ham,  1952)  and  the  supersonic 
area  rule  of  Hayes  (1947)  and  Lomax  (1955).  The  supersonic 
area  rule  shows  that  the  pressure  disturbance  for  a  complex 
three-dimensional  configuration  can  be  reduced  to  the  calculation 
of  the  pressure  disturbance  due  to  an  equivalent  body  of  revolution, 
provided  that  the  position  of  interest  is  sufficiently  far  from 
the  body.  Whitham’s  theory,  on  the  other  hand,  describes  an 
asymptotic  flow  behavior  at  a  distance  sufficiently  far  from  a 
body  of  revolution.  This  asymptotic  flow  obeys  the  geometric 
acoustic  laws.  That  is,  the  flow  disturbances  are  linearly  pro¬ 
portional  to  a  local  F  function,  which  is  related  to  the  shape 
of  the  body  and  to  the  flow  conditions.  The  values  of  the  F 
function  are  constant  along  the  characteristic  curves  emitted 
from  the  body.  Consequently,  the  sonic  boom  signature  at  a 
distance  far  from  an  aircraft  can  be  calculated  from  the  F  function 
of  an  equivalent  body  of  revolution  of  the  aircraft. 

In  experimental  investigations  of  the  sonic  boom  in  wind 

tunnels,  it  is  usually  necessary  to  use  very  small  models  in 

order  to  obtain  direct  measurements  of  the  far-field  pressure 

signature  in  the  vicinity  of  a  wind  tunnel  wall  (see,  for  example, 

2.1 


Carlson  (1964)).  With  these  small  models,  inaccuracies  with 
respect  to  model  contours,  vibration  of  model,  boundary  layer 
development,  interferences  of  the  sting  supports,  small  non¬ 
uniformities  of  the  free  stream,  etc.  usually  arise.  The  present 
proposed  wind  tunnel  testing  method  is  directed  towards  an 
alleviation  of  this  very  small  model  restriction. 

The  present  proposed  method  is  based  on  large  models  in 
wind  tunnels,  where  only  the  near  or  the  mid-field  is  simulated. 

By  measuring  the  pressure  distribution  at  the  vicinity  of  the 
wind  tunnel  wail,  it  is  possible  to  determine  the  signature  of 
the  sonic  boom  at  large  distances  in  the  outer  mid-field  or  in 
the  far-field.  This  method  is  expected  to  have  great  significance 
in  making  wind  tunnel  tests  more  feif able  by  avoiding  the  use 
of  extremely  small  models  is  usually  done  today. 

It  is  well  known  that,  in  the  near  field  of  an  aircraft, 
the  flow  field  is  fully  three-dimensional.  The  supersonic  area 
rule  may  not  be  applied  there  and,  hence,  Whitham’ s  supersonic 
projectile  theory  may  not  be  employed  directly.  Moreover,  the 
flow  disturbances  in  the  near  field  may  not  be  generally  * 
described  by  Whitham* s  asymptotic  relations  even  within  the 
assumption  of  the  linearized  supersonic  flow  (Pan  1970  asb, 
1971).  Indeed,  recent  wind  tunnel  experiments  (Morris,  Lamb  and 
Carlson,  1970)  have  "shows  that  Whitham* s  theory  does  not  give 
good  predictions  6f  pressure  signature  shape  in  the  near  field 
especially  at  large  Mach  numbers  and  at  large  angles  of  attack. 


Consequently,  the  calculation  of  flow  disturbances  in  a  not-so- 
far  field  from  an  aircraft  or  the  extrapolation  of  a  sonic  boom 
pressure  signature  from  a  known  pressure  disturbance  at  a  nearer 
field  must  be  examined. 

In  this  section,  based  on  the  well-known  quasi-linear 
assumptions,  we  shall  extend  the  current  sonic  boom  theory 
to  the  near  field .  We  shall  be  particularly  concerned  with 
the  non-axisymmetric  and  non-geometric-acoustic  effects  in  the 
flow  field.  Based  on  the  new  extended  theory,  we  shall  present  a 
new  wind  tunnel  testing  method  for  sonic  boom  based  on  large 
models. 

2.1.2  EXTENSION  OF  THE  CURRENT  SONIC  BOOM  THEORY 

As  discussed  in  the  preceding  sub-section,  the  current 
sonic  boom  theory  is  limited  to  predict  the  pressure  dis¬ 
turbance  at  a  distance  sufficiently  far  from  a  supersonic  air¬ 
craft.  In  this  sub-section,  we  shall  develop  a  theory  to  be 
applicable  almost  in  the  entire  flow  field. 

GENERAL  CONSIDERATION 

Let  us  consider  a  steady,  homogeneous,  irrotational, 
supersonic  flow  over  a  three-dimensional  configuration  such  as 
an  aircraft  (Figure  2.1.1).  The  flow  disturbance  generated  by 
the  aircraft  in  the  supersonic  flow  propagates  along  character¬ 
istic  surfaces  to  the  far  field  in  the  downstream.  It  is  clear 
physically  that  the  flow  disturbance  in  the  far  field  due  to  the 
presence  of  this  aircraft  is  equivalent  to  the  flow  disturbance 
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generated  from  a  streamtube  of  upstream  radius  R,  or  any  other 
streamtubes,  enclosing  this  aircraft  in  the  nearer  field. 

A  streamtube  of  upstream  circular  cross  section  with 
radius  R  is  chosen  arbitrarily  but  not  too  close  to  the  air¬ 
craft  (Figure  2.1.1),  the  flow  disturbance  over  this  stream- 
tube  is  generally  weak  and  this  direamtube  is  of  quasi-cylindrical 
shape  with  a  mean  radius  R.  The  flow  field  over  this  streamtube 
can  be  treated  by  the  well  known  linearized  supersonic  flow 
theory.  The  shape  of  this  streamtube  and  the  flow  disturbance 
over  this  streamtube  at  R  can  be  represented  by  a  superposition 
of  various  multipole  distributions.  By  defining  a  Whitham 
type  F  function  in  terms  of  these  multipole  distributions,  we 
can  obtain  the  F  function  for  this  streamtube  of  upstream 

radius  R,  which,  in  turn,  can  be  related  to  the  various  flow 

* 

disturbances  at  R. 

The  propagation  of  the  known  F  function  representing 
flow  disturbance  from  one  streamtube  to  another  one  further 
afield  may  be  treated  basted  on  Whitham*  s  hypothesis  on  the  im¬ 
provement  of  characteristics;  that  is,  the  values  of  the  F  function, 
are  constant  along  bicharacteristics  in  the  three-dimensional 
flow.  Thus  the  new  F  function  is  obtained  for  a  streamtube 
further  afield.  However,  the  new  F  function  may  geneially 
have  multiple  values  in  certain  regions -of  its  arguments ; 
these  are  due  to  the  intersection  of  characteristic  surfaces  in 
the  physical  space  where  values  of  physical  quantities  cease 
to  be  unique.  This  failure  of  the  linearized  supersonic  flow 
theory  as  a  description  of  the  flow  is  known  to  be  remedied  by 
the  presence  of  shock  surfaces.  The  positions  of  shock  surfaces 


may  be  determined  by  the  usual  simple  geometric  property;  that 
is,  to  the  first  order  of  shock  strength,  the  shock  surface 
bisects  the  angle  between  the  two  intersecting  characteristic 
surfaces-  Having  fixed  the  positions  of  the  shocks,  the  new 
F  function  becomes  single  valued,  which  actually  is  the 
Whitham  type  local  F  function  of  the  corresponding  stream- 
tube  further  afield.  Finally,  from  the  new  single  valued  F 
function  for  this  streamtube  further  afield,  the  corresponding 
flow  disturbances  may  be  obtained. 

GOVERNING  EQUATION  AND  ITS  SOLUTION 

Let  us  choose  the  body  axis  to  be  the  x-axis  coinciding 

with. the  free  stream  direction  (Figure  2.1.2).  Enclosing  this 

body  we  choose  a  co-axial  quasi-cylindrical  streamtube  with  an 

* 

arbitrary  upstream  radius  R  (R=0  corresponds  to  a  general  point¬ 
nosed  slender  body).'  On  the  surface  of  this  streamtube  the  flow 
is  assumed  to  be  disturbed  at  x  =  0  where  the  origin  Of  the  x-axis 
is  located.  The  r-axis  perpendicular  to  the  x-axis  is  the 
radial  coordinate  and  0  is  the  polar  coordinate  measured  counter¬ 
clockwise.  from  the  vertical  downward  r-direction.  We  further 
assume  that  the  entire  problem  is  symmetric  with  respect  to  9  =  0 
(or  v)  plane. 

) 

Let  the  free  streets  velocity  and  Mach  number  be  U  and 

? 

M  respectively,  and  at  a  general  point  (x,r,8)  the  local 
velocity  be  (U+ Uu,Uv,Uw)  ]  The  flow  is;  assumed  to  be  irrotational, 
hence  the  velocity  disturbances  u,v,  andw  may  be  deduced  from 
a  velocity  potential  On  the  linearized  flow  theory,  4>  satisfies 
the  equation, 

2.6 


aratubeand  Coordinate  System 


(2.1.1) 


d2$  I  d$  ,  1  d2$  2  d$2 

— _  + - +  ■ — — .  -  «  — -pc-  = 

dr  r  dr  r2  d 9Z  bxZ 


0, 


where  p  =  (M2-  l)i/2 


A  general  solution  of  Eq. (2.1.1)  which  represents  a 
disturbance  propagating  downstream  from  the  quasi-cylinder  is, 
(Ward  1955), 


4> 


r- 

=  -  \  cos  t9  I  -S_££ 
n  -pR 


[(x-t)2-p2r2 
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(2.1.2) 


disturbance  velocity  components  deduced  from  Eq.  (2.1.2)  are. 


u  =  - 

/>  AC 

x-pr 

r 

nn  1 

L 

n 

tub 

ny  j 

-pR 

[(s-t)2  -  ' 

COS 

ne 

x-pr 

r  x-t 

v  = 

-pa  r 

-  sV]W 

and 

n 

x  'pr 

ne  j 

J-pR 

W  s: 

sin 

Here, 

hn<*> 

=  cosh  j 

n  cosh~*tJi 

5 

(2.1.3) 


(21.4) 


(2.1.5) 


(2.1.6) 


and  fn(t)  are  the  multipole  distributions  and  may  be  related 
to  the  shape  of  the  quasi-cy lindr ica 1  st reamt ube  (Lomax  and 


He  as  let  ,  1956)*.  We  shall  solve*  for  f (t )  in  terms  of  the  flow 


disturbances  later.  The  pressure  disturbance  relates  to  the 
velocity  disturbances  by  the  linearized  Bernoulli's  equation 


(Ward  1955), 


P-P. 


00 


(2.1.7) 


F  FUNCTION  AND  ITS  FOURIER  COMPONENTS 


By  changing  the  integration  variables  and  by  replacing 
x-gr  =  y  -  pR  in  Eqs.  (2.1.4)-(2.1.6),  we  have 


—  I 


cos  n  0 


1 


hn(1+^T)  fn(*-(3R)dt 

I7T 


JCy-t)  (y-t  +  2pr>y 


(2.1.8) 


cos  n  0 


y-t  +  er 
r 


«r> 

[(y-t)  (y-t  +  2pr)]  1/2 


(2.1.9) 


w 


I 


sin  n  0 


i  V1*  dt 

r  [(y-t)  (y-t +  2pr)] 1/14 


(2.1.10) 


Here,  y  is  the  charact eristic  paraaeter  of  a  linearized  character 
istic  curve  from  the  surface  of  the  streaatube  of  upstream 
radius  R. 

To  define  a  Whitham  type  F  function,  we  make  a  far-field 
approximation,  i.e.,  gr/y  »  1  and  write  velocity  components  in 
Fourier  series  expansions.  The  expressions  (2.1.8)  to  (2.1.10) 
are  reduced  respectively  to, 


u  cos  n  Q 
n 


(2.1.11) 


u  = 


L 


n 


(231- ) 


172  l_,  Fn(y)  cosn 


e.. 


n 


v  =  ^  vn  cos  n  0  =  — — 3, 


(23D 


i?2  Fn(y) 


cos n 0  (2.1.12) 


n 


=  1 


w  =  ^  w  sin  n  0  = 

n 


r(23r) 


2/2  n  Gn(y)  sinne  (2.1.13) 


Here, 


y  f ' (t  -  3R)  dt 


?«(y)  ■/  -^^72 


(2.1.14) 


being  the  Fourier  components  of  the  F  function  of  the  stream- 

-  1 

tube  with  upstream  radius  R 


F(y,0) 


FnXy)  cos  n  0 


(2.1.15) 


and 


o 


Fn(t)  dt. 


(2.1.16) 


It  may  be  noted  that  for  bodies  of  revolution  (n=o) ,  the 

above  relations  are  those asymptotic  linear  relations  of  Whitham 

(1352);  for  bodies  of  revolution  at  small  angles  of  attack 

(n=o  and  1),  the  above  relations  reduce  to  those  of  Siegelaan 

(196/7).  It  may  also  be  pointed  out  that  in  the  above  relations 

—1/2  -9/2 

u  and  v  vary  as  r  '  and  w  varies  as  r  to  the  far  field; 
hen.de,  in  the  far-field,  w  may  be  neglected  in  comparison  with 


u  or  v-  Consequently,  u,  v  and  P  may  be  related  only  to  the 
F  function  of  an  equivalent  body  of  revolution.  This  is  the 
basis  on  which  the  current  sonic  boom  calculations  are  made. 

RELATIONS  BETWEEN  THE  F  FUNCTION  AND  FLOW  DISTURBANCES 

It  has  been  shown  by  Pan  (1970  a,b)  that  it  is  possible, 
within  the  linearized  supersonic  flow  theory,  to  obtain  exact 
relations  between  the  flow  disturbances  and  the  corresponding 
local  F  function  on  an  arbitrary  streamtube  enclosing  an  axisymmetric 
body.  These  relations  are  valid  in  the  entire  flow  field  of  an 
axisymmetric  slender  body  in  the  supersonic  flow.  Now  we  shall 
obtain  similar  relations  for  the  flow  over  a  non-axisymmetrie 
quasi-cylindrical  streamtube  of  upstream  radius  R  (Figure  2.1.2). 

At  large  R,  these  relations  reduce  to  Eqs.  (2.1.11)-(2.1.13) . 

As  stated  previously,  the  shape  of  the  streamtube  of  up¬ 
stream  radius  R  is  represented  by  the  superposition  of  the 
various  multipole  distributions  fE,  where  fn  may  be  detertoihed 
in  terms  of  the  flow  disturbances  on  the  surface  of  the  Strdam- 
tube  by  the  boundary  conditions  of  Eq.  (2.1.1).  By  setting  r=R 
in  Eqs. (2.1.3)- (2.1.5)  [or  Eqs.  (2.1.8)-(2.1.10)j ,  using- 
Eq. (2.1.7),  and  expanding  the  flow  disturbances  in  Fourier 
series  in  6  [See,  Eqs.  (2.1.11)-(2.1.13)] ,  the  multipole  dis¬ 
tributions  may  be  related  to  the  Fourier  components  of  the 
flow  disturbances  by  the  following  integral  equations;  at  (x,R), 

n  n  J  £  <x-t)  (x-t  +  20E)]1/i! 


J 


(2.1.17) 


(2.1.18) 


J 


X-t+BRhn(W+1)  fn.(t-PR>  dt} 
R  £x-t ) (x-v  + 20R)] X/^ 


and 

r  n 

W  —  |  ""  "*  . r  1  l  /o 

n  JQ  R  [  (x~t)  (x-t  +  2$R)]1/j£ 


(2.1.19) 


These  are  integral  equations  for  f^  of  fQ  and  can  be  solved. 

Either  f  or  f’  can  be  found  in  terms  of  the  Fourier  components 
n  n 

of  the  flow  disturbance  P  .u  ,v„  or  w  .  Solutions  have  been 

n  n  n  n 

obtained  and  reported  by  Pan  (1971).. 

By  substituting  the  solutions  of  fn  and  into  the  de¬ 
finitions  of  Fn  and  Gn  [Eqs.  (2.1.14)  and  (2.1.16)],  Fn  and  Gfl 
can  be  written  in  teams  of  the  flow  disturbances  (Pan,  1971),  at 
R 
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=  < 


Pn(x) 


>  - 


PR 


A 
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pnit) 


-^(t) 


>  Sa(W>  dt' 
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-  '  i  I 


V«  Tn(^r)  dt> 


and 


tiGa(x) 


x 

+  j 

bR  JL 


V«  sn(W>  *• 


(21.20) 


(2.1.21) 


(2.1.22) 


1 7$  =  wn(x> 

Here  S  and  T  are  resolvent  Kernels  of  K*  and  j',  respectively. 

n  H  '  XI  « 


ana  are  the-  differentiations  of  KQ  and  Jn,  respectively 
with  respect  to  their  arguments,  where 

*/2  2 

03/2  r  h  (x  cos  z-fl)dz 


u  „  (2  +  x  cos  z) 


Jn(x)  = 


(x  cos  z  +  1)  K  (x  cos  z  +  1)  dz 


(2  +  x  cos  z) 


(2.1.24) 


Sn,Tn’Kn  and  Jnf  have  been  computed  on  the  UTSI  computer  IBM 

1130  and  are  tabulated  respectively  in  Tables  2.1.1  to  2.1.4 

for  s  from  0  to  2.0  and  for  n  from  0  to  5.  TS  .T^K*  and  J* 

1  o  o  o  o 

were  denoted  respectively,  as  S,  S^,  K* ,  and  previously  by 
Pan  (1970a)]. 

The  Fourier  components  of  the  flow  disturbances  may  also  be 
expressed  in  terms  of  the  corresponding  Fn  or  Gn,  at  R  (Pan  1971) 


-Un(x 


L—i 


(2pR) 


x 

[f  (x )  +  —  f  F  (t)K'(S^)  dtl  (2.1.25) 

L  PR  Jn  a  pR  J 


vn(x)  - -  1) 

(2pR)' 


Fn(x)  + 


--  j  F  (t)  j'{— )  dt]  (2.1.26) 

^  "  11  AD  ! 


wn(x)  = - 

R(2pR)A/ 


G  (x)  +—  J  G  (t)  k!(—  )  dt  (2.1.27) 
.  pR  o  n  n  pR  J 


Since  the  streamtube  of  upstream  radius  R  was  chosen 


TABLE  2.1.1  FUNCTIONS  Sn(x) 
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Values  of  data  beyond  this  point  cannot  be  accurately  obtained  at  present. 
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arbitrarily,  the  relationship  between  the  flow  disturbances  and 


Fn  or  Gn  [Eqs.  (2. 1. 20)- (2 . 1 . 22)  or  Eqs.  (2 . 1.25)- (2. 1.27) ]  are 
valid  for  any  R  with  in  the  linearized  supersonic  flow  theory. 

As  R-*»,  these  relations  reduce  to  the  asymptotic  linear  relations, 
Eqs. (2. 1. 11)- (2. 1. 13) ,  which  for  n=0  correspond  to  Whitham's 
asymptotic  relation  (1952).  The  second  terms  in  Eqs.  (2.1.20)- 
(2.1.22)  and  Eqs.  (2. 1.25)- (2. 1.27)  represent  the  corrections 
of  the  asymptotic  relations  in  the  nearer  field. 


PROPAGATION  OF  THE  DISTURBANCE 


Since  the  relations,  Eqs.  (2. 1.25)- (2. 1.27) ,  aie  valid 
where  the  flow  disturbances  are  weak,  we  may  now  write  the 
expressions  for  flow  disturbances  at  any  distance  r  in  terms 
of  corresponding  local  F  function  at  r. 


1  A£ 

Poo 

-u(y) 


r  }  ,  _  I 

— 72"  F(y,0  )  +  —  cos  n  0Q  j  **n(t)  V^)dt 
3r)  'I  B**  n  )  pr  J 


Br  J 
(2.1.28) 


v  (y)  = 


- g  TTo  F(y  ,9)+—  Y*  cos  n  e  f 

(2pr)1/2  -  °  Br  7*  °  "o 


o  J 

o  8r 


w(y)  = 


r  (2gr) 


(2.x. 29) 


-  r\  Tn~  [G(y,9)  +  ~  ^  cosne  J  G  (t) 

6*r/2  |deL  ^  O 


/  G  (t)K' (^)dt 
o  pr 


Based  on  Whitham’s  hypothesis  on  the  improvement  of  the 
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linear  theory  (1952) ,  the  flow  disturbance  in  steady  supersonic 
two-dimensional  flow  propagates  along  an  improved  two-dimensional 
characteristic  which  is  at  an  angle  with  respect  to  the  flow 
direction  equal  to  the  Mach  angle.  In  three-dimensional  flow, 
the  two-dimensional  characteristic  is  replaced  by  a  bi-characteristic 
direction  and  a  characteristic  surface  orientation.  A  bi-characteristic 
direction  is  one  which  is  at  an  angle  with  respect  to  the  flow  dir¬ 
ection  equal  to  the  Mach  angle.  A  characteristic  surface  orienta¬ 
tion  is  one  for  which  the  angle  between  the  normal  and  the  flow 
direction  is  the  complement  of  the  Mach  angle.  A  bicharacteristic 
line  is  an  integral  curve  for  the  bicharacteristic  direction.  A 
characteristic  surface  is  a  surface,  the  orientation  of  which  is 
everywhere  the  characteristic  orientation.  A  general  bicharacteristic 
line  through  a  given  point  is  far  from  unique;  however,  a  characteristic 
surface  through  a  given  smooth  line  in  space  is  locally  unique.  Any 
envelope  of  characteristic  surface  is  also  a  characteristic  surface. 

To  a  given  bicharacteristic  direction  at  a  point  corresponds  a 
unique  characteristic  orientation,  for  which  the  bicharacteristic 
direction  is  tangent  to  the  characteristic  surface.  Through  a 
given  non-singular  point  in  a  given  characteristic  surface  passes  • 
a  unique  bicharacteristic  line  in  a  characteristic  surface  (Hayes 
and  Probstein,  1966). 

For  the  present  problem  of  the  linearized  supersonic  ilow,  the 
characteristic  surfaces  are  surfaces  of  revolution  with  respect  to 
the  axis  (Ferri,  1954).  Through  a  general  point  (x,r,9)  on  a 
characteristic  surface,  the  bicharacteristic  direction  may  be  found 
from  the  tangential  direction  to  the  characteristic  surface  and  to 


the  local  Mach  cone. 


,  idM- S,,'V«  ••• 


2  2  1/2 

Ax  _  (q  -v  )  '  (l-f-u)cosu,  -  (1-Ht)v  sin  it 
Ar  (q2  -  v2)1^2  v  cos  n  +  (q2  -  v  )siny 


2  2  1/2 

rAS  „  (q  -  v  )  cos  u  -  v  sin  u 
=  w 


Ar 


(q  -  v  )  v  cos  y  +  (q*  -  v sin  y 


with  y  and  q  being  the  local  Mach  angle  and  the  local  non- 
dimensional  velocity  respectively. 

By  using  the  definitions  of  y  and  q,  and  neglecting  the 
second  and  higher  order  terms  of  the  flow  disturbances,  we  have 
the  differential  equations  of  the  bicharacteristics 


dx 

Sr 

d0 

dr 


=  p  +  u  -  M2  (v  +  pu) 


1$ 


m 

r 


(2.1.31) 


(2.1.32) 


To  obtain  a  set  of  parametric  equations  of  the  bicharacteristics, 
we  may  substitute  Eqs.  (2.1.28)-(2.1.30)  into  Eqs. (2. 1.31)  and 
(2.1.32)  and  perform  the  integrations  on  the  line  y=  constant 
and  eQ  =  constant , 


r 

x  =  3(r- R) +y-k  F(y,0o)(ri/2- R1/2)  -  /  L(y,0o; 


r)dr , 


(2f.  1.33) 


and 

r 

0  =  0Q  +  M(y,eo;r)dr,  (2.1.34) 

R 

with 

k  =  2"1/2(7+l)  M4p~3/2  (2.1.35) 

y 

Uy,e0;r)  =  VCOsne  f  f  (t)K’(^)dt 

°  (3(2pr)3/2  ~n  JQ  pr 

+  M2 (v  +•  pu)  (2.1.36) 

and 

M(y,0o;r)  =  pw/r .  (2.1.37) 

Here,  y  and  9Q  are  the  coordinates  x  and  0,  respectively,  of  a 
point  on  the  cylindrical  surface  R.  y  =  constant  and  0Q  =  constant 
define  a  bicharacteristic  curve  from  the  said  point;  on  this 
bicharacteristic  curve,  F(y,0Q)  =  constant.  For  an  axisymmetric 
flow  field  (n=0), there  is  no  0-component  velocity  and  0  remains 
constant;  Eqs.  (2.1.31)  and  (2.1.32)  reduce  to  Whitham’s 
characteristic  differential  equation,  and  'q. (2.1.33)  reduces 
to  the  improved  characteristic  equation  obtained  by  Pan  (1970  a,b). 

After  knowing  F(y,0Q)  at  R,  we  may  integrate  Eqs.  (2.1.33) 
and  (2.1.34)  to  obtain  F(x,r,0)  at  a  cylindrical  surface  further 
afield  at  r.  The  value  of  F  at  a  particular  point  (x,0)  at  r 
is  equal  to  the  value  of  F(y,0Q)  at  a  point  (y,0Q>  which  is  on 
the  same  bicharacteristic  curve  as  the  point  (x,0)  at  r.  Due 
to  the  ©-component  disturbance  the  bicharacteristic  curve  from 


(y,e  )  is  generally  not  lying  on  the  same  &  plane  (see, 

Figure  2.1  3a).  The  bicharacteristics  which  pass  through  the 
straight  line  0  =  constant  on  the  surface  r  originate  on  the 
surface  at  R  from  the  curve  0Q  =  0Q(y).  Consequently,  to 
determine  F(x,r,0)  on  a  line  0  =  constant  at  r,  we  have  to 
determine  the  curve  $Q  =  9  ( y)  at  R.  On  the  plane  of  symmetry 
0=0  (or  it),  where  w  =  0,  all  bicharacteristic  curves  remain  in 
the  same  plane;  hence  F(x,r,0)  is  determined  only  by  F(y,0)  at 
R  (see ,  Figure  3.1.3b). 

The  shape  of  the  F  function  obtained  at  r  are  generally 
distorted  and  different  from  the  original  F  function  at  R. 

The  F  function  on  a  0  =  constant  line  at  r  is  generally  a 
multi-valued  function  of  its  argument  x.  The  multi-valued 
regions  are  due  to  the  intersection  of  characteristic  surfaces 
in  the  physical  space  where  values  of  physical  quantities 
cease  to  be  unique.  This  failure  of  the  linear  theory  as  a 
description  of  the  flow  is  known  to  be  remedied  by  the  presence 
of  shock  surfaces  (Courant  and  Friedrichs,  1948). 

It  is  well  known  (Courant  and  Friedrichs,  1948)  that  the 
shocks,  to  a  first  order  in  strength,  can  be  determined  by  a 
simple  geometric  property*,  that  is  a  shock  surface  bisects  t*“~ 
angle  between  two  intersecting  characteristic  surfaces.  As 
shown  by  Whitham  (1552)  the  F  function,  at  least  in  the  far 
field,  gives  a  rough  description  of  the  flow  pattern,  since  it 
shows  whether  the  characteristics  are  converging  in  compression, 
where  a  shock  will  appear,  or  diverging  in  expansion.  The  shock 
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position  at  any  distance  r  may  be  determined  from  the  F  function  j 

?, 

at  R.  For  an  axisymmetric  flow,  Whitham  (1952)  obtained  a  | 

15 

relation  called  the  "area-balance-rule"  which  states  that  the  | 

lobes  cut-off  on  each  side  of  the  F  function  by  a  straight 

segment  which  determines  a  shock  position  must  be  equal  in 

area.  The  slope  of  the  straight  segment  depends  on  the  distance 

r.  In  the  following,  we  shall  follow  Whitham* s  procedure  to 

obtain  general  relations  for  determining  shock  positions. 

Suppose  that  a  general  shock  intersects  9  =  constant  at  r 

! 

at  the  point , 

» 

I 

x  =  p(r-R)  -  G (r , 0 )  (2.1.38) 


and  the  bicharacceristics,  Ea.  (2.1.33),  specified  by  y^  and  y2 
(y2  >  y-j^)  on  0o  =  0o(y)  at  R  intersect  the  shock  at  this  point 
(see,  points  c  and  d  in  Figure  2.1.3).  The  bisection  of  the 
angle  between  the  characteristic  surfaces  by  the  shock  surface 
requires  that 


2  &  =  I  k  r1/2 
3r  2 


F(yl’eo(yl))+  F(y2’Vy2}) 


L(yl'eo(yl)  *r)  +  L(y2,eo(y2)  ,r)j 


(2.1 .39) 


On  the  other  hand,  elimination  of  x-s(i*~R)  from  the  equations  of 
the  shock,  Eqs.  (2.1.38)  and  of  the  bicharacteristic,  Eq. (2.1.33), 
give 


•naMugqj)^^ 


and 


r 

G(r,0)  =  kF(y1.6?o(y1)).r1/2-R1/2j  +  J 

J  R 


L^r0o(yi);r)dr"yr 

(2.1.40) 


G  (r ,  0 )  =  k  F(y2’e0^y2^ 


rV2  .  r1/2 


']  +  f  .L(y2^0(y2);r)dr“y2' 

J  JR 

(2.1.41) 


G  and  r  as  functions  of  y^  or  y2  can  be  solved  from  Eqs.  (2.1.39)- 
(2.1.41).  and  then  the  relation  between  y^  and  y2  is  obtained. 

The  relation  between  r  and  y^  (and/or  y2)  is  obtained  by 
eliminating  G  from  Eqs  (2.1.40)  and  (2.1.41), 


F(y2>0o(y2^  "  F(yi,0o(yl)  * 


t 


y2  '  yl 


-L 

f 


=  k(r1/2-H^2)  + 


L(y2,e0(y 2);r)  -L(y1,eo(y1);r)j 
P(y2,fl0(y2»-  F(y1,eo(y1)) 


dr 


-1 


(2.1.42) 


For  larger,  Eq.  (2.1.42)  may  be  simplified  to  a  relation  obtained 
by  Whithaa  (1952)  for  R  =  0, 


F(y2’go(y2)}  ~  F(yl>9o(yl)) 
y2  “  yl 


(2.1.43) 


The  geometric  interpretation  of  Eq.(2.1.43)  is  that  the  slope 
of  the  straight  segment  cd  joining  the  points  y^  and  y2  of  the 
F  function  curve  at  R  relates  to  r  only.  For  an  arbitrary  r, 
the  slope  of  the  segment  cd  relates  not  only  to  r  but  to  y^»y2 
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and  F.  For  different  shocks,  the  slopes  of  the  segments  are 
different . 


To  determine  the  relation  between  y^  and  or  the  position 
of  cd  on  the  F  function  curve  at  R,  the  well  known  "area- 
balance-  rule"  , 

y<> 


J  F(yieo(y)  )dy  =  2  (y2~y2)  F(yl,0o(yl))  f  F(y2,0o(y2^]  5 


(2.1.44) 


is  found  to  be  replaced  by  the  following  relation, 


/  f(y  eo(y))dy  =  A  <y2‘yi) 


[F(y1,80(y1))  -  F(yz,eo(y2): 


r 


F(yr£1o(yl))!  J  fL(yl>6,o(yl);r)  "  L(y2,eo(y2);r) 

j  R  L 


dr 


r 

r 


[  F(y,9  ly)Wy  —  j  L(y,0  (y);r, 

v.  '**J* 


dr 


(2.1.45) 


Now  Eq. (2.1.45)  together  with  Eq. (2.1.42)  determines  the 
positions  of  y^  and  y2  on  F(y10o(y)}  for  fixed  r.  These  relat  '^ns 
may  be  used  to  determine  the  positions  of  other  shocks.  For 
example,  the  position  of  the  front  shock  ab  (Figure  2.1.3)  may 
be  determined  by  setting  F(yi’0o(yi>>  -  0  at  che  point  a  on 
the  F  function  curve.  We  may  note  here  that  Eqs.  (2.1.42)  and 
(2.1.45)  are  so  complex  that  no  explicit  solution  for  r(y1.y2» 


2.26 


. . . 


possible.  Results  may  only  be  obtained  by  a  numerical  iterative 
procedure  by  using  Eqs.  (2.1.43)  and  (2.1.44)  as  a  first 
approximat ion . 

After  having  determined  the  positions  of  shock  waves,  the 
local  F  function  at  r  becomes  a  single  valued  function  except 
with  a  finite  number  of  discontinuities.  Then  the  corresponding 

flow  disturbances  may  be  determined  from  the  local  F  function 

* 

and  its  Fouriex*  components  by  Eqs.  (2.1. 28)-(2. 1.30)  or  the 
relations  (2. 1.25)- (2. 1.27)  by  replacing  R  by  r. 

2.1.3  WIND  TUNNEL  TESTING  METHOD  BASED  ON  LARGE  MODELS 

As  discussed  in  the  subsection  2.1.1,  in  experimental  in¬ 
vestigations  of  sonic  boom  in  wind  tunnels  it  is  usually 
necessary  to  use  very  small  models  in  order  to  obtain  direct 
measurements  of  the  far  field  pressure  signature  in  the 
vicinity  of  a  wind  tunnel  wall.  The  present  new  method  is 
based  on  large  models  in  wind  tunnels,  where  only  the  near  or 
the  mid-field  is  simulated.  By  measuring  the  pressure  dis¬ 
tribution  at  the  vicinity  of  wind  tunnel  wall  it  is  possible 
to  determine  the  signature  of  sonic  boom  at  large  distances  in 
the  far-field.  In  this  subsection,,  we  shall  present  a  theoretica 
study  of  the  new  wind  tunnel  testing  method. 

GENERAL  CONSIDERATION 

A  three-dimensional  model  is  tested  in  a  circular  cylindrica 
supersonic  wind  tunnel  with  radius  R  (Figure  2.1.4a).  The  flow 
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is  assumed  to  be  steady,  homogeneous,  irrotational  and  inviscid. 
It  is  well  known  that  the  disturbance  emitted  from  the  model 
is  reflected  at  the  wind  tunnel  wall,  in  order  that  the  stream¬ 
lines  adjacent  to  the  wall  may  be  always  tangent  to  the  wall. 
Hence,  the  three-dimensional  reflected  pressure  disturbance  can 
be  measured  at  the  wind  tunnel  wall.  If  the  flow  field  over  the 
test  model  inside  the  wind  tunnel  can  be  described  by  the 
linearized  supersonic  flow  theory  and  if  the  cumulative  non-' 
linear  effects  can.be  neglected  within  the  wind  tunnel,  the 
flow  field  near  the  wall  may  also  be  described  by  the  linearized 
supersonic  flow  equations.  It  is  possible  to  relate  the  incident 
disturbance  with  the  reflected  disturbance  at  the  wall  by 
specifying  proper  boundary  conditions  on  the  wall  (Pan,  1970c, 
1971). 

In  free  flight,  however,  streamlines  over  the  model  may 
be  distorted  freely  according  to  the  emitted  disturbance.  Any 
circular  cross-sectional  streamfube  in  the  free  stream  may  be 
distorted  into  a  quasi-circular  cylindrical  shape.  The  quantity 
of  flow  in  the  wind  tunnel  of  radius  R  is  equal  to  the  quantity 
of  flow  confined  in  a  streamtube  of  the  same  upstream  radius  R 
enclosing  this  model  in  free  flight  (see,  Figure  2.1.4b). 

Based  on  the  linearized  supersonic  flow  theory,  the  actual 
flow  disturbance  at  the  distance  R  in  free  flight  equals  to 
the  incident  disturbance  on  the  wind  tunnel  wall  with  the 
same  radius  R. 

After  knowing  the  actual  pressure  disturbance  at  a  stream- 
tube'  with  upstrei'm  radius  R  enclosing  the  model,  we  may  now 
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apply  the  extended  theory  developed  in  the  preceding  subsection 
to  extrapolate  this  known  disturbance  to  any  distance  in  further 

<  j 

afield. 

WIND  TUNNEL  WALL  REFLECTION 

In  this  subsection,  we  shall  obtain  the  relation  between 

3 

*  j 

the  actual  free  flight (incident)  pressure  disturbance  and  the  j 

J 

measured  (reflected)  pressure  disturbance  both  on  the  cylindrical 

wind  tunnel  wall.  We  assume  that  the  flow  field  may  be  described 

by  the  linearized  supersonic  flow  theory.  The  pressure  (Mach) 

waves  generated  by  the  model  propagated  outward  and  downstream  ; 

I 

and  are  reflected  at  the  wind  tunnel  wall;  the  reflected  wave 

i 

I 

will  not  interact  with  the  model.  Unlike  the  reflection  of  a 

i 

plane  wave  at  a  rigid  plane  surface  which  gives  a  constant  re-  j 

flection  factor  2.0,  the  reflection  of  a  curved  wave  from  a 
cylindrical  rigid  surface  may  not  be  described  by  a  simple  con¬ 
stant  reflection  factor.  Since  the  flow  field  behind  a  curved 
wave  is  generally  not  uniform,  the  flow  field  behind  a  curved 
reflected  wave  depends  on  both  the  incident  curved  wave  and  the 
distance  from  the  wave  front. 

Referring  to  Figure  2.1.4a  with  the  coordinates  described 
in  Figure  2.1.2,  the  disturbance  velocity  potential  $  in  the 
wind  tunnel  satisfies  Eq. (2.1.1).  Equation  (2.1.1)  together 
with  the  uniform  free  stream  boundary  conditions  may  be  solved 
by  a  technique  of  Laplace  transform,  i.e.. 


jr 
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(2.1.46) 


5(s,r,0)  =  —  I  4>(x,r,0)e  SX//^Rdx. 

3R  o 

The  general  solution  to  the  transformed  governing  equation 
Eq. (2.1.1)  is  given  by  (Ward,  1955). 

5(s,r ,0)  =  V  cos  n  0  Fa  K  (— )  +  B  I  (—  >1,  (2.1.47) 

4*  L  n  n  R  nnP.J 

where  K  and  I  are  modified  Bessel  functions,  and  A  and  B 
n  n  n  n 

are  functions  of  s  only.  [Readers  should  not  be  confused  by  the 
modified  Bessel  functions  used  here  with  the  functions  defined 
by  the  Eqs  (2.1.23)]. 

In  Eq. (2.1.47)  the  terms  in  Kn  which  represent  waves  travel¬ 
ing  outward  from  the  model  and  downstream  are  the  incident  waves 
with  respect  to  the  wall,  while  the  terms  in  In  which  represent 
waves  traveling  in  both  directions  are  related  to  the  waves  re¬ 
flected  from  the  wall.  The  relation  between  An  and  BQ  may  be 
determined  by  the  boundary  condition  on  the  wall,  (i.e., 
d$/dr  =  0  at  r  =  R) , 

B  =  -  AK'(s)/l'(s).  (2.1.48) 

n  n  n  n 

The  function  A  may  be  determined  from  the  boundary  conditions 
n 

on  the  model  surface',  we,  however,  shall  leave  it  as  an  arbitrary 
function. 

Now  the  general  solution  in  (s,r,0)  reduced  to 


PR  =  [cosne 


(2.1.51; 


n 


with  PRn  being  the  Fourier  components  of  PR,  and  from  Eqs. 
(2.1.49)  and  (2.1.50), 


~  Ws> 


1  - 


Kn(s)  Vs> 
i;(s)  ^(s) 


(2.1.52) 


On  the  otner  band,  it  is  easy  to  show  that  the  Fourier  components 
of  the  free  flight  (incident)  disturbance  pres sore  on  the  wall  is 


PT  (s,R)  =  -  —  A  K  (s) 
ln  PR  n  n 


(2,1.53) 


Hence,  we  may  relate  the  measured  pressure  with  the  free  flight 
pressure  on  the  wind  tunnel  wall, 
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i  I  rfTwfa 


1 


PR  (s,R) 

Kn 


=  PT  (s,R) 
in 


1 


K>>  In(s) 
i;<S)  Kn(s) 


(2.1.54) 


We  may  also  express  the  free  flight  pressure  in  terms  of  the 
measured  pressure  and  use  the  well,  known  Bessel  function 
relations  (Watson,  1966), 


PT  (s,R) 
xn 


-  Pv  (s ,  R)  s  K  (s) 
2  Rn  n 


(s) 


+  In-l<E)j 


(2.1.55) 


Now  the  relations  between  PT  (x,R)  and  P0  (x,R)  may  be 

in  «n 

obtained  by  performing  the  inverse  Laplace  transform, 

PT  (x ,  R)  =  -  ,  PR  (x,R) 
xn  2  [ 

=  0 


(2.1.56) 


Here,  Rn(x)  the  reflection  functions,  are  the  inverse  Laplace 
transform  of  the  Rn(s), 


Rn(s)  =  1  -  s  Kn(s) 


fws> 


4-  I 


n-1 


(2.1.57) 


Since  the  laplace  inversing  of  Rn(s)  is  not  known  at  the  present, 
exact  values  of  R^te)  for  all  n  and  x  cannot  be  obtained. 

However,  for  most  practical  cases,  the  argument  of  Rn(x)  is 
usually  small,  and  we  may  obtain  an  asymptotic  expression  of 
Rn(x)  for  small  x  by  inverting  the  asymptotic  expression  of 
Rn(s)  for  large  s.  The  result  is 


2.33 


2.34 


of  the  extrapolation  method  and  list  the  relevant  formulas  in 
each  procedure. 

A  three-dimensional  pressure  distribution,  P^(x,0.R)  is 
measured  a--  a  circular  cylindrical  wind  tunnel  of  radius  R  (for 
example,  by  using  a  scanning  system  suggested  by  R.  C.  Bauer 
of  ARO,  Inc.  to  this  author,  1971).  This  three-dimensional 
signature  is  expanded  in  Fourier  series  in  6, 


PR(x,0,R)  =  ^  P^fe.R)  cosne; 
n 


the  Fourier  coefficients  PRn  of  P^  can  be  obtained. 

The  corresponding  Fourier  coefficients  PT  (x,R)  of  the 

xn 

incident  pressure  disturbance  Pj(x,0,R)  and  Pj(x,0,R)  itself 
can  be  found,  from  Eq. (2.1.56), 


PT  (x,R)  =-  [pR  (x, R)  -  —  f  PR  (t , R)  R(— )  dt 
In  2  I.  pr  J  pR 


and  from  Eq. (2.1.59), 


Pj (x, R,0)  =  ^  Pj^(x, R)  cosn0. 
n 


Pj(x,R,0)  is  a  three-dimensional  pressure  disturbance  on  a 
streamtube  with  upstream  radius  R,  which  corresponds  to  the- 
wind  tunnel. 

By  using  the  pressure  disturbance  Pj(x,R,0)  or  Pj  (x,R), 
we  can  calculate  the  local  F  function  of  this  streamtube;  from 
Eq.  (2.1.20), 


I 
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1 


Fn(y,R) 


=  (20R) 


1/2 


y 

p  (y )  ~  f 

.  n  PR  J 


ln 


(t)  Sn(^i)  dt 
n  PR  -» 


*  ^ 


and  from  Eq.  (2.1.15), 


F(y,e0,R)  =  Fn(y,R)  cosneQ. 
n 

F(y,eQ,R),  representing  the  disturbance  generated  by  the  stream- 
tube  with  upstream  radius  R  or  by  the  tested  model  in  free 
flight  at  R,  propagates  to  further  afield  r  along  bicharacter¬ 
istics  given  by  Eq.  (2.1.33), 

r 

x  =  p(r-R)  +  y  -  k  F(y,0o)  (r^-R1^2)  -  J  L(y,eQ;r)dr 

R 

and  Eq.  (2.1.34), 
r 

9  =  eQ  +  J  M(y,0o;r)dr. 

R 

After  shock  positions  were  fixed  on  the  F  function  by 
Eqs.  (2.1.42)  and  (2.1.45),  the  new  single  valued  F(x,r,0) 
is  obtained  at  r.  This  new  F  function  is  then  expanded  in 
Fourier  series. 


F (x , r , 9 )  =  ^  Fn (x , r )  cos  n  9 . 
n 

Corresponding  to  the  Fourier  components,  Fn(x,r),  the  Fourier  com¬ 
ponents  of  the  pressure  disturbance,  Pn(s,r),  are  found  from 
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. . . — . . . 


Eq.  (2.1.25). 


pn  (x  > r  )  = 


(2f3r) 


1/2 


Fn(x,r)  + 


if 

erJ 


Fn(t,r)  Kae 


s±>  dtl 

3r  j 


Finally,  the  three-dimensional  pressure  disturbance  P(x,r,0) 
at  r  is  obtained  by  summing  up  all  Fourier  components, 


P(x,r,e)  = 


n 


PR  (x , r )  cos  n  9 . 


2.1.4  EXAMPLES  AND  DISCUSSIONS 

$ 

Several  numerical  calculations  of  near-field  pressure 
signatures  were  performed  to  compare  the  results  based  on  the 
present  analysis  with  those  based  on  the  current  sonic  boom  cal¬ 
culation  and  with  some  available  wind  tunnel  measurements.  Typical 
examples  are  presented  and  discussed  in  the  following. 

Figure  2.1.6  shows  a  comparison  of  the  shifts  of  character¬ 
istics  at  different  r  for  different  points  x  on  the  axis  of  a 
6.46°  half-angle  cone-cylinder  body.  The  length,  t,  of  the  cone 
portion  of  this  body  is  taken  to  be  1.0  and  the  free  stream 
Mach  number  is  1.41.  In  this  figure,  yw  is  the  shift  of  the 
asymptotic  characteristic  curve  based  on  Whitham's  theory  (1952) 
from  the  linear  characteristic  curve;  y^,  is  the  shift  of  the 
characteristic  curve  based  on  the  present  analysis  (Pan,  1970a). 


*This  author,  Dr.  Y.  S.  Pan,  wishes  to  acknowledge  Misters  K.  T. 
Wang  and  M.  O.  Varner,  Research  Assistants,  for  their  valuable 
assistance  in  the  numerical  computations. 
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VIGURE  2.1.6  Cone-cylinder  body  and  shifts  of  character-  ' 

istics  at  different  r  ! 

I 
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For  this  example,  the  present  characteristics  are  located 
generally  ahead  of  the  asymptotic  characteristics.  At  r  =  0.5, 
the  differences  of  the  shifts  are  as  large  as  20  percent  of 
Whitham' s  values.  As  r  increases,  the  percentage  of  the  dif¬ 
ference  decreases . 

Figures  2.1.7a  -  2.1.7b  show  the  near-field  pressure 
signatures  at  r/t  =0.5,  1.0  1.5,  2.0,  2.5,  3.0,  4.0  and  5.0 

from  the  axis  of  the  cone-cylinder  body  described  in  the  preceding 
section.  Comparisons  show  the  differences  between  the  signatures 
of  the  present  analysis  and  of  the  current  sonic  boom  theory. 

The  differences are  due  to  the  present  correction  of  the  linear 
asymptotic  relations  between  the  flow  disturbances  and  the  local 
F  function  and  due  to  the  shift  of  the  characteristic  curves. 

The  peak  pressures  predicted  by  the  present  analysis  are  generally 
lower  than,  and  located  ahead  of,  those  predicted  by  the  current 
sonic  boom  theory.  This  prediction  is  qualva*. ‘xely  consistent 
with  some  of  the  near-field  experimental  Oi  -  r/ations  (for  example, 
Morris,  iAmb  and  Carlson,  1970). 

In  Figures  2.1.7b  -  2.1.7h,  numerical  pressure  data  are 
also  presented.  These  data  were  obtained  by  Kutler  (1971) 
using  a  shock-capturing  finite-difference  approach  (Kutler  a^d 
Lomax,  1971).  The  present  signatures  are  generally  in  excellent 
agreement  with  the  numerical  data  except  at  the  neighborhood  of 


^Numerical  data  generously  provided  by  Dr.  Paul  Kutler  of  NASA 
Ames  Research  Center  are  gratefully  acknowledged. 
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the  shock  waves.  The  difference  near  the  shock  waves  is  due  to 
the  fact  that,  in  the  shock-capturing  finite-difference  approach, 
the  shock  waves  which  occur  are  spread  over  a  few  mesh  intervals. 
Therefore  the  exact  positions  of  the  shock  waves  cannot  be  exactly 
predicted  by  this  numerical  method. 

Figure  2.1.8  shows  the  shifts  of  characteristics  at  r  =  1.0 
from  the  same  cone-cylinder  body  at  Mach  numbers  1.2,  1.41, 

2.0  and  3.0.  Generally,  the  shift  of  characteristics  increases 
with  increasing  Macn  number  except  at  the  cone  portion.  In 
this  region,  the  F  function  is  positive  and  the  shift  decreases 
with  increasing  large  Mach  number.  Figure  2.1.9a  -  2.1.9d  show 
the  comparison  of  pressure  signatures,  based  on  two  different 
analyses,  at  r/t  =  1.0  of  the  same  cone-cylinder  body  at  several 
different  Mach  numbers.  Because  of  the  larger  shifts  of  character¬ 
istics  at  larger  Mach  numbers,  larger  shifts  of  pressure  signatures 
appear  at  larger  Mach  numbers.  The  difference  of  the  pressure 
strengths  decreases  as  the  Mach  number  increases. 

Figure  2.1.10  and  2.1.11  demonstrate  the  extrapolation 
of  an  axisymmetric  pressure  signature.  From  Kutler's  numerical 
data  (1971)  at  r/t  =  1.5,  a  corresponding  F  function  is  obtained 
and  is  shown  in  Figure  2.1.10  (the  dotted  data  represents  the 
F  function  based  on  Whitham’s  asymptotic  relation).  This  F  function 
is  propagated  to  r/t  =  2,0,  2.5,  3.0,  4.0,  and  5.0.  The  corres¬ 
ponding  pressure  signatures  at  different  r/t  are  obtained  from 
different  loca?  F  functions  and  ai’e  shown  in  Figures  2.1.11a- 
2.1, He.  Comparisons  of  the  extrapolated  pressure  signatures 
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with  Kutler's  numerical  data  are  again  in  good  agreement  except 
at  the  neighborhoods  of  the  shock  waves. 

Figure  2.1.12)  shows  the  pressure  signatures  of  a  slender 
body  of  revolution  at  r/-t=4.2  and  M  =  1.2.  Experimental  data 
(Barger,  1968)  are  compared  with  the  results  of  the  two  different 
theoretical  calculations.  The  difference  between  the  two 
theoretical  curves  is  obvious.  By  observation,  one  cannot  tell 
which  curve  is  in  better  agreement  with  the  experimental  data, 
because  the  streamwise  positions  of  the  experimental  data  with 
respect  to  the  body  are  usually  not  specified.  Based  on  Whithara’s 
theory,  Barger  (1968)  developed  a  procedure  to  design  a  body  of 
revolution  from  a  specified  signature.  Pan  and  Varner  (1971) 
performed  a  similar  calculation  by  using  the  present  analysis. 
Examples  indicate  the  different  body  radii  obtained  from  the  two 
different  analyses.  The  differences  are  most  apparent  at  the 
rearward  portion  of  the  bodies. 

Figure  2.1.13-2.1.18  demonstrate  an  extrapolation  of  a 
measured  three-dimensional  pressure  signature.  A  three-dimensional 
measured  pressure  signature  is  assumed  in  a  form,  PR  =  P^  +  PRj  cos  0 
and  is  shown  Figure  2.1.13.  A  free  stream  Mach  number  of  M  =2.0 
and  a  wind  tunnel  wall  radius  of  R  =  0.5  are  assumed.  The 
magnitude  of  the  assumed  pressure  disturbance  has  been  exaggerated 
to  apply  the  usual  linearized  approaches.  Following  the  pro¬ 
cedures  described  in  Subjection  2.1.3,  the  corresponding  incident 
(free  flight)  pressure  disturbance  is  calculated  (n  =  0  and  1)  and 
is  shown  in  Figure  2.1.14. 
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Present  Analysis  (Pan,  1970a) 

Current  Sonic  Boom  Theory  (Whitham,  1952) 
Experiments  (Barger,  1968) 


FIGURE  2.1.12  Calculated  and  measured  pressure  signatures 


From  the  incident  pressure  disturbance  of  its  Fourier  com¬ 
ponents  (n  -  0,1)  the  Fourier  components  of  the  corresponding 
F  function  are  first  calculated.  By  summing  up  the  Fourier 
components  (n  =  0,1),  the  local  F  function  is  obtained  and  is 

shown  in  Figure  2. 1. 15.  Since  the  kernel  functions  S  and  S,  in 

o  1 

the  formula  for  calculating  Fn  and  Fn  are  of  different  signs  but 
with  the  same  order  of  magnitude,  the  near-field  effects  of  the 
F  function  are  almost  negligible  in  this  example.  For  more  com¬ 
plicated  signatures  where  more  Fourier  component  terms  are  re¬ 
quired,  the  near-field  effects  on  the  function  may  be  important. 

The  disturbance  represented  by  the  F  function  is  extra¬ 
polated  to  the  downstream  at  r  =  10.0.  The  local  single  valued 
F  function  at  r  =  10.0  is  obtained  by  following  the  procedures 
described  in  Subsection  2.1.3  and  is  shown  in  Figure  2.1.16. 

The  F  function  at  r  =  10.0  is  then  expanded  in  a  six-term  Fourier 
series  in  0(n  =  0  to  5).  This  six-verm  Fourier  series  is  generally 
a  good  representation  of  the  f  function  except  at  certain  positions 
(e.g.  shock  waves)  where  the  F  function  is  discontinuous  in  0. 

From  each  Fourier  component  of  the  F  function,  six  Fourier  com¬ 
ponents  of  the  corresponding  pressure  disturbance  are  obtained. 

By  summing  up  the  six  Fourier  components,  the  corresponding 
three-dimensional  pressure  signature  is  obtained  and  is  shown  in 
Figure  2.1.17.  In  this  figure,  the  pressure  signatures  near  the 
leading  shock  wave  are  fixed  approximately  from  the  six-term 
Fourier  representation  of  the  pressure  signatures  which,  for 
0  =0°,  90°  and  180®,  are  shown  by  dashed  lines  in  the  same  figure. 
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FIGURE  2,1.18  Comparisons  of  pressure  signatures 


If  we  take  more  terms  in  the  Fourier  series  in  the  F  function, 
and  consequently  in  the  pressure  disturbance,  the  actual  signature 
will  approach  the  signature  indicated  by  the  solid  lines. 

In  Figure  2.1.18,  the  pressure  signatures  at  6  =  0°,  45°, 

90°,  135°  and  180°  obtained  previously  are  compared  with  those 
obtained  directly  from  the  current  sonic  boom  theory.  In  this 
example,  the  current  sonic  boom  theory  over-estimates  the  peak 
pressures  for  0°  0  £  90°  as  much  as  20  percent  but  under¬ 

estimates  the  peak  pressures  for  90°  <  0  ^  180°  as  much  as  50 
percent.  Although  no  general  quantitative  cone lu5,ion  on  the 
importance  of  the  near-field  effects  can  be  obtained  from  this 
single  example,  this  example  does  show  that  the  near-field  ef¬ 
fects  may  be  important  for  extrapolating  complicated  near-field 
three-dimensional  signatures. 

2.1.5  CONCLUSIONS 

The  current  sonic  boom  theory,  which  is  based  on  the  well 
known  "Supersonic  Area  Rule"  and  Whitham’s  "supersonic  projectile 
theory" ,  has  been  extended  to  take  into  account  the  three- 
dimensional  near-field  effects.  The  present  extended  theory 
permits  us  to  calculate  the  flow  disturbances  in  the  almost 
entire  flow  field  of  a  steady  supersonic  flow  over  a  three- 
dimensional  configuration. 

Based  on  the  present  theory,  a  new  wind  tunnel  testing 
method  based  on  large  models  is  developed.  This  new  method  per¬ 
mits  us  to  use  the  models  of  usual  sixe  in  wind  tunnels  to  deter¬ 
mine  the  sonic  boom  on  the  ground  by  measuring  the  pressure 
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signature  at  tne  wind  tunnel  wall.  This  method  has  great 
significance  in  making  wind  tunnel  tests  more  reliable  by  avoid¬ 
ing  the  use  of  extremely  small  models  as  is  usually  done  today. 

2.1.6  RECOMMENDAT IONS 

THEORETICAL  WORKS 

It  is  recommended  to  compute  a  more  complete  set  of  the 

/  / 

kernel  functions  S„,T  .K  and  J  which  are  required  for  the 

a  n  n  n 

application  of  the  present  theory  to  the  practical  sonic  boom 
calculations  and  to  the  extrapolation  of  pressure  signatures 
measured  in  the  wind  tunnels.  Only  part  of  these  kernel 
functions  have  been  computed  and  presented  in  this  report. 

It  is  recommended,  based  on  the  present  theory,  to  develop 
a  FORTRAN  computer  program.  More  three-dimensional  calculations 
are  required  to  compare  with  experimental  data  and/or  with  the 
calculations  based  on  the  current  sonic  boom  theory. 

It  is  recommended  to  develop  a  sonic  boom  theory  valid  at 
large  Mach  numbers.  the  range  of  applicability  of  the 

present  quasi-linear  approach  is  quite  limited,  a  theory  which 
takes  into  account  non-linear  effects  at  large  Mach  numbers 
should  be  developed  and  should  be  of  great  importance  in  the  near  future. 

EXPERIMENTAL  WORKS 

It  is  recommended,  based  on  the  present  method,  to  develop 
wind  tunnel  test  procedures  to  produce  suitable  three-dimensional 
pressure  signatures  for  extrapolation  to  the  far  field  that  is 
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2.2.  FUNDAMENTAL  THEORIES,  APPLICABILITY  AND  EXTENSIONS 

2.2.1.  The  Near- Field  Flow  Pattern  of  an  Inclined  Slender 

Body  of  Revolution 

by  Dr.  S.  N.  Chaudhuri,  and  LTSI  student  Sarat  Praharaj 

Summary 

In  the  present  paper  Whitham’s  far-field  theory  of 
supersonic  flow  pattern  has  been  generalized  for  the  near¬ 
field  points  of  an  inclined  body  of  revolution.  Using  the 
correct  expre.^-  on  for  the  characteristics  (as  far  as  the 
first-order  theory  is  concerned)  it  is  shown  that  the  new 
F-function  which  describes  the  near-field  flow  pattern  is 
dependent  on  the  following  besides  the  body  geometry:  the 
distance  from  the  axis  of  the  body,  the  Mach  number,  the 
angle  of  attack  and  the  azimuthal  plane  angle.  We  have 
obtained  closed  form  results  for  the  new  F-function  for  any 
smooth  body  of  revolution.  The  recurrence  formulas  given  are 
convenient  for  computer  programming.  The  pressure  signatures 
for  near-field  positions  have  also  been  calculated.  They 
can,  of  course,  no  longer  be  related  to  the  new  F-functions 
by  simple  formulas.  The  routine  method  given  to  calculate 
the  fund  mental  F-function  and  the  pressure  signature  at  near¬ 
field  p^j-nts  of  non-axisymmetric.  bodies  will  be  useful  in  the 
preliminary  design  study  of  equivalent  bodies  of  revolution. 


2.64 


LIST  OF  SYMBOLS 


=  Angle  of  attack 
=  ,(M2  -  i)V2 

=  Velocity  potential 
=  Perturbation  velocity  potential 

=  Axial-flow  perturbation  velocity  potential  (First  Order) 

=  Cross-flow  perturbation  velocity  potential  (First  Order) 

=  Whitham's  function  for  the  points,  whore 
pr/x-pr  is  large 

Ft  y,r ,M,R(y) ,a,0}  =  Whitham’s  function  at  any  point 


a 

3 

$ 

P 

G 


O 

CPl 

F(y) 


f(t) 

g(t) 

M 

Po 

AP 


(x) 

(r) 

(9) 


Ii(y) 

U 


Strength  of  source  (singularity)  distribution 
Strength  of  doublet  (singularity)  distribution 
Free  stream  Mach  number 
Free  stream  static  pressure 

Pressure  in  excess  of  undisturbed  (free  stream) 
static  pressure 
total  velocity 

total  axial  velocity  component 

total  radial  velocity  component 

total  azimuthal  velocity  component 

Distance  of  a  point  on  the  characteristic  from  the 

axis  of  the  body 

Radius  of  the  body  of  revolution  at  y 
Free  stream  speed 

Distance  of  any  point  measured  from  the  nose 
Distance  of  any  point  from  the  nose,  where  the 
characteristic  produced  meets  the  axis. 
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Introduction 


For  an  axisyrametrical  flow  past  a  slender  body  of 
revolution  Whitham  (1952)  developed  an  elegant  modified  theory 
which  essentially  retains  the  simplicity  of  the  linearized 
theory  as  given  by  Lighthill  (1952)  and  others  but  remedies  the 
failure  of  the  linear  theory  as  a  description  of  the  flow 
pattern.  The  fundamental  hypothesis  of  Whitham  is  that  the 
linearized  theory  gives  valid  first  approximation  to  the 
flow  everywhere  provided  that  in  it  the  approximate  character¬ 
istics  are  replaced  L.  sufficiently  good  approximation  to 
the  exact  ones.  This  hypothesis  has  been  amply  substantiated 
by  checks  detailed  in  Whitham*  s  paper  (Whitham,  1952'). 

The  following  are  some  of  the  important  assumptions 
made  in  Whitham* s  theory  (Whitham,  1952): 

(i)  The  body  is  slender  and  pointed  at  the  nose,  with 
the  front  shock  attached.  (It  may  be  remarked 
here  that  even  if  these  conditions  are  not  satisfied, 
Whitham* s  theory  can  still  be  used  to  deduce  the 
behaviour  of  the  flow  at  large  distances  from  the 
axis  of  the  body). 

(ii)  Whitham* s  discussions  of  the  flow  pattern  is  limited 
strictly  to  the  behaviour  at  large  distances  from 
the  axis  of  the  body.  In  particular,  the  Whitham 
function  F(y)  which  is  fundamental  to  the  whole 
theory  and  is  the  most  important  function  associated 
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i 

I 


with  flow  past  a  body  of  revolution  is  obtained 

only  for  far-field  prsitions.  The  condition  is 

also  valid  for  front  shocks  at  any  distance.  This 

simplifies  the  mathematical  analysis  considerably 

and  makes  the  F-function  dependent  only  on  y,  the 

distance  of  any  point  from  the  nose  where  the 

characteristic  produced  meets  the  axis. 

(iii)  Apparently  a  third  restrictive  condition  in  Whitham's 

theory  is  that  of  axial  symmetry.  However,  Ward (1949) 

+ 

has  shown'  that  the  flow  becomes  axisymmetrical 
when  (3r/(x-pr)  is  large  (See  Fig.  2.2J}  Here  p=  (M2-l)ly/2  j 

with  M  as  the  free-stream  Mach  number  and  r  is  the  ] 

distance  of  a  point  on  the  characteristic  from  the 
axis  of  the  body.  The  quantity  (x-pr)  is  the 
linearized  from  of  the  characteristic  variable  and 
measures  the  distance  from  the  nose  at  which  the 
characteristic  starts.  This  Ward's  condition 
is  clearly  satisfied  at  large  distances,  but  it 
is  also  true  at  points  on  the  front  shocks  because 
the  appropriate  characteristic  surfaces  arise  so 
very  close  to  the  nose.  Hence  the  results  for  the 
front  shock  and  all  the  Wbitham  theory  at  large 
distances  apply  unchanged  t.  ?  r.on-axisymraetrical 
slender  body.  From  thos  •  Wnixham  concludes  (Whitham,  1952) 
that  it  is  reasonable  to  expect  that  his  results 

+  This  also  folxows  from  the  result  established  in  the  present 
paper  as  will  be  shown  later. 
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for  large  distances  apply  to  the  supersonic  flow 
past  any  finite  body.  Mathematically  this  leads 
to  the  interesting  result  that  even  for  non- 
axisymmetrical  bodies  the  simplified  Whitham  F- 
function  can  be  used  for  far-field  points. 

We  shall  now  examine  the  above  assumptions  critically 
with  a  view  to  its  application  to  the  calculation  of  near¬ 
field  flow  pattern  of  an  inclined  slender  body  of  revolution. 

The  first  assumption  is  essential  to  our  theory  because 
the  first-order  potential  equations  are  used.  Strictly  we 
do  not  use  the  entire  linearized  equation  in  our  analysis 
as  it  would  be  had  the  x-axis  been  aligned  with  the  free- 
stream  direction.  This  was  first  noted  by  Lighthill  (1948) 

The  consequences  of  this  distinction  has  been  considered  in 
detail  by  Van  Dyke  (1952)  and  we  will  discuss  it  later. 

The  second  assumption  is,  in  general,  obviously  in¬ 
valid  for  near-fields.  It  is,  however,  still  true  for  the 
front  shock  system  originating  very  close  to  the  nose  but 
not  for  other  points  of  the  body.  Whitham  (1952)  replaces  the 
approximate  straight  characteristics  x  -  f3r  by  a  better 
approximation,  y(x,r)  given  by  equation  (10)  in  (Whitham,  19."2). 

Whitham  justifiably  simplifies  the  "extremely  complicated" 
expression  (10)  for  his  far-field  theory  and  obtains  the 


Incident ly,  there  is  an  error  in  equation  (10).  The  second 
integral  on  the  right-hand-side  of  the  equation  should  be 
multiplied  by  (-m2)  and  not  by  (-2M2) .  This,  however,  does 
not  affect  Whitham’ s  results  as  he  has  not  used  equation 
(10)  in  his  paper. 
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well-known  results  for  the  equation  of  the  characteristics 
[equations  (12)  and  (13)  of  Whitham’s  paper  (1952)) 


x  =  (3r  —  k  F(y )  r  '  +  y 


where 

„  (y-t) 


k  =  2-1/2(7  +  l)  M4  p-3/2 


Whithara  notes,  however,  that  expression  (10)  of  his  paper  (Whit ham, 1952) 
still  provides  the  correct  approximation  to  the  characteristic 
curves  near  the  body  and  should  be  used  for  near-and  mid- 
field  calculations. 

Using  the  correct  expression  for  y(x,r)  as  given  in 

equation  (10)  of  the  above  reference  it  will  be  observed  that  the 

new  F-function  will  be  dependent  on  y,r,M  and  R(y)  (see  list 

of  symbols)  for  axisymmetrical  flows.  We  have  denoted  this 

function  i>y  F  {  y  ,r,M,R(y)}  .  For  inclined  flows  there  are 
s 

two  more  variables  on  which  F-function  will  depend,  namely, 
the  angle  of  attack  a  and  the  azimuthal  plane  angle  a. 

We  have  obtained  closed  form  results  for  the  new  F- 
function  ..or  any  smooth  body  of  revolution.  The  recurrence 
fox*mulas  given  are  convenient  to  use  for  computer  programming. 


A  O  I 

'In  Whitham's  paper  a  -  (M  -1)  , which  we  have  replaced  here 

by  3  in  order  to  use  the  Greek  letter  a  for  angle  of  attack. 
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It  is  shown  that  if  the  shape  of  the  meridian  curve  is  ex-  ■ 
pressible  as  a  polynomial,  all  the  integrals  are  related  by 
simple  recurrence  formulas  and  ultimately  depend  on  one 
standard  form.  It  is  a  fortunate  fact  that  any  function, 
which  is  continuous  in  a  closed  interval,  can  be  uniformly 
approximated  within  any  prescribed  tolerance,  over  that  interval, 
by  some  polynomial.  This  follows  from  the  well-known  theorem, 
which  states  that  for  a  continuous  function  y(x)  in  an  inter¬ 
val  (a,b)  there  will  be  a  polynomial  f(x)  such  that  |f(x)-y(x)| 

^  e  in  (a,b)  for  an  arbitrary  positive  e. 

Once  our  polynomial  f(x)  has  been  determined  so  that 
it  satisfactorily  approximates  the  given  meridian  curve  R(x) 
over  a  certain  interval  (a,b)  it  is  easy  to  see  that  integration, 
being  essentially  a  smoothing  process,  will  involve  lesser 
error  than  the  derivations  in  f(x)  and  R(x)  (Hildebrand,  1956).  It  is  of 
importance  to  emphasize  here  the  fact  that  the  integrands 
representing  the  F-functions  are  of  the  type  G(x)  f(x), 
where  G(x)  are  given  functions  and  f{x)  is  the  approximated 
polynomial.  Therefore  the  closed  form  solutions  obtained 
will  introduce  even  smaller  error  than  by  using  mechanical 
quadratures. 

The  third  restrictive  condition  in  Whitham’s  theory  is 
that  of  axial  symmetry.  As  mentioned  earlier  the  far-field 
flow  pattern  is  not  significantly  affected  by  removing  this 
assumption  of  axial  symmetry.  But  is  is  obvious  that  we 
have  to  consider  the  effects  of  cross-flow  for  bodies  at  an 


angle  of  attack  in  the  calculation  of  the  flow  pattern  at 
near-field  points.  This  has  been  done  by  a  distribution  of 
doublets  with  their  axes  perpendicular  to  the  body  center- 

line.  The  strength  of  the  doublets  per  unit  length  is  found 
to  be  proportional  to  the  cross  section^area  if  we  replace  the 
exact  tangency  condition  by  an  approximate  one.  This 
approximation  is  exactly  similar  to  the  first  order  axial 
flow  pioblem. 

To  illustrate  the  type  of  F-curve  that  results  from  the 
present  theory  we  have  calculated  the  new  function  F  for 
two  bodies.  Their  shapes  are  given  by 

(a)  R(x)  =  0.1  x(l-x),  0  £  x  £  1 

=  0  ,  x  >  1 

(b)  R(x)  =  0. 1  {  1-  (1-x)3} ,  0  ^  x  ^  1 

=  0.1,  x  >  1 

A  comparison  of  the  near-  and  mid  field  F-function  with  the 
far  field  F- functions  bring  out  the  following  important 
facts: 

(i)  As  expected  F-function  very  nearly  coincides  with 
the  F-function  for  points  close  to  the  nose. 

(ii)  At  other  stations  there  are  significant  differences 
with  changes  in  r,  especially  the  peaks  and  troughs 
of  F-function. 
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(iii)  For  axisyrametric  flows  F-f unction  is  a  weak  function 
of  M  upto  about  3.  This  again  is  to  be  expected 
from  the  slender  body  theory, 
liv)  For  inclined  flows  the  effect  of  Mach  number  is 
more  pronounced. 


improved  First-Order  Theory  for  Supersonic  Flow  Past 
Inclined  Body  of  Revolution 


Consider  uniform  supersonic  flow  past  a  body  of  re¬ 
volution  inclined  at  an  angle  of  attack  a  [see  Fig.  2.2.1J. 

The  boundary  condition  at  the  body  is  simplified  by  choosing 
a  cylindrical  coordinate  system  (x,r,0)  aligned  with  the  bcdy 
axis.  We  resolve  the  tree-stream  velocity  into  axial  and 
cross-flow  components,  as  shown  in  Fig.  2.2.2.  The  shape  on  the 
body  is  defined  by  its  continuous  meridian  curve  c  -  R(x) . 

For  moderate  supersonic  Mach  numbers  the  flow  is  assumed 
to  be  isej .tropic  and  irrotational  throughout  so  that  there 
exists  a  velocity  potential  <j>(x,r ,0).  The  equation  of  motion 
is  "then 


<*Z  -  #  «»,  +  <“2  -  «?>  "rr  *  <•* ' -  \ ^  \ ♦«.  +  <»2  ;  hr 


2  H  r^0^'r0  ~  2  Gx  r^xr  0 

r  r 


(2.2. la) 


where  the  speed  of  sound  a  is  related  to  its  free  stream 


value  a_  by 


a  =  ao  "  *  (*x  'r  4r  +  "2  *e  -  U  > 

z  r 


(2.2. lbN 


The  suffixes  have  been  used  to  denote  partial  differentiations. 
Introducing  a  perturbation  potential  U^(x,r,e)  we  write 


$(x,r,e)  =  U[x  cos  a  +  r  sin  a  cos  6  +  0(x,r,0)]  (2.2.2a) 


The  total  axial,  radial  and  azimuthal  velocity  components 
are  then  given  by 


q(x)  -  *x  -  u(cosa  +  **>> 
q^  =  4>r  =  U(sin  a  cos  0  +  ffr). 


,  (2.2.2b) 


q(0)  =  ~  ^9  =  u(“  sin  a  sin  9+  - 
r  r 

Substitutingeq.  (2.2.2a)  into  equations  (2.2. 1)  gives  the  exact  equation 
in  the  perturbation  velocity  potential.  We  need  not  write 
the  full  equation  for  cur  purpose.  However,  in  order  to 
bring  out  the  distinction  between  first-order  and  linearized 
equation  for  this  inclined  flow  problem  we  reproduce  the 
perturbation  equation  partially  retaining  only  terms  linear 
in  tf : 


*rr  +  1  K  +  4  -  -  ‘‘xx  sin2cl  +  *rr  sin2=t  cos2s 


11  2  2 

+  (-  — cL)sin  a  sin  9  +  sin  2  a  cose 

_  uv  „  xr 

r  r 

-  —  rf  „  sin2a  sin  6  — —  <rf..sin2a  sin29 


r  —  &  sin  ct  sin  6  H — o  sin^ct  sin2e 

„  r  a  (i 

r  r 

•r  the  nonlinear  terms  in  0  (2.2.3a) 

To  this  we  add  the  following  boundary  conditions: 

(i)  All  perturbations  should  vanish  upstream  of  the 
body  so  that 

0  (0,1,0)  =  0(Q,r,G)  =0  (2.2.3b) 

Xv 

(ii)  The  flow  s‘r  uld  be  tangential  to  the  surface  of 
the  body  so  that 

0  (x,R,0)  +  sina  cost)  =  R'[cosa  +  d  (x,R,0)],  (2.2.3c) 

X  XL 


where 

R'  _  dR(x) 
dx 

Actually,  the  first  upstream  conditions  should  be  supple¬ 
mented  by  the  suock  condition  along  the  bow  wave.  However, 
this  refinement  is  not  necessary  for  the  order  of  accuracy 
sought  here. 

The  First-Older  Problem  and  Its  Justification 

The  complete  perturbation  equation  (2.2.3a)  including  the 
non-linear  terms  is  equivalent  to  the  origin"- 1  non-linear 
equation  of  motion,  Eqs. (2.2.1).  Consequently,  simplifying 
assumptions  must  be  made  in  order  to  solve  it.  The  well  known 
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theory  of  Karman  and  Moore  (1932 )  or  axially  symmetric  flow  and 
of  Tsien  (1938)  or  inclined  flow  is  based  upon  the  assumption 
that  the  entire  right-hand  side  of  Eq.  (2.2.3a)  can  be  neglected 
leaving  the  wave  equation.  Its  solution  will  be  termed  the 
first-order  potential  (Van  Dyke ,  1951 )  and  is  denoted  by  <p  (x ,  r ,  e ) . 

It  is  clear  from  Eq.  (2.2.3a)  that  the  first-order  and 
linearized  problems  are  the  same  where  the  free  stream  is 
aligned  to  the  x-axis,  i.e.,  where  a  =  0.  But,  as  mentioned 
earlier,  Lighthill  (1948)  first  noted  that  wbaia  /  0  the  linearized 
problem  is  not  the  same  as  the  first-order  solution  because 
we  are  neglecting  the  linear  terms  in  cp  in  Eq.  (2.2.3a).  The 
justification  advanced  for  this  approximation  for  inclined 
flows  follows  from  the  slender-body  theory  in  which  the  fol¬ 
lowing  are  some  of  the  important  assumptions  generally  made: 

(i)  approximate  pressure  relation;  (ii)  approximate  tangency 
condition;  and  (iii)  small  angle  of  attack.  Hence,  Lighthill (1948) 
and  Laitone  (1947)  conclude  that  the  exact  solution  of  the  first 
order  equation  is  fruitless  and  that  first-order  theory  is 
incapable  of  yielding  results  more  exact  than  those’  of  slender 
body  theory. 

The  first-order  problem  is,  therefore, 

<Prr  +  +  ~  ^2®xx  =  0  (2.2.4a) 

q>(0,r,{?)  =  v-ox(0,r,9)  =  0  (2.2.4b) 
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q»  (x,R,0)  +  sina  cos 0  R  [ cosa  +  m  (x ,  R ,  f.) )  I . 

J.  A 


(2.2.4c) 


The  first-order  problem  is  satisfied  exactly  by  setting 


tp(x,r,0)  =  q>Q(x,r)  cosa  +  cpjCx.r)  sina  cose  (2.2.5) 


The  first  term  corresponds  to  the  axial  component  of  free 
stream  velocity  and  the  second,  to  the  cross-flow  component 
(Fig.  2.2.2).  We  shall  assume  that  a  is  small  so  that  cosa^  1 
and  sin  a  ~  a. 

The  i irst-order  problem  is  thus  separated  into  two 
completely  independent  problems.  For  the  axial  flow, 


orr 


f* 


0 


(2.2.6a) 


;0(°,r)  =  <pox(o,r)  =  0  (2.2.6b) 

mor(x’R)  =  r/[1  +<Pox(x,R)J  (2.2.6c) 


and  for  the  cross-flow. 


!'lrr 


+ 


(<r>1r/r)  “ 


(c^/r2)  - 


f32cp 


XX 


=  0 


(2.2.7a) 


cp^Ojr)  =  =  0  (,2.2.7b) 

1  +  q>lr(x,R)  =  R,<plx(x,R)  (2.2.7c) 


The  general  solutions  of  Eqs.  (2.2.6a)  and  (2.2.7a)  satisfying  up¬ 
stream  conditions  (2.2.6b)  and  (2.2.7b)  are  given  in  (Von  Karman, 
1932)  and  (Tsien,  1938). 
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(2.2.8a) 


x-pr 

©o(x,r)  =  -  J  X(t,f3r)  f(t)  dt , 
o 

x-pr 

on,  (x,r)  =  —  /  (x-t)  X(t,pr)g(t)  dt,  (2.2.8b) 

o 

where 

-1/2 

X(t,pr)  =  { (x-t)2  -  p2r2}  (2.2.8c) 

♦ 

The  corresponding  first  derivatives  are 

x-Pr 

<Pox  =  -  f  X<t,pr)  f'(t)  dt,  (2.2.9a) 

o 

x-pr 

cp  =  -  /  (X-t)  X(t ,f3r)  f'(t)  dt,  (2.2.9b) 

r 

o 

x-pr 

tp,„  =  —  /  (x-t)  X(t  ,pr)  g'  (t)  dt,  (2.2.9c) 

o 

x-pr 

cplr  = - ^  /  (x-t)2  X(t,6r)  g'  (t)  dt  (2.2. 9d) 

6r  Q 

The  functions  f  and  g  are  strictly  determined  from  the 
tangency  conditions,  Eqs.  (2.2.6c)  and  (2.2.7c)  which  leads  to 
Volterra  integral  equations,  which  can,  in  general,  only  be 
solved  numerically.  For  our  order  of  accuracy  as  used  by 
Whithara  (1952)  we  may,  neglect  the  term^  compared  with  cos  a 
in  the  exact  tangency  condition  of  Eq.(2.2.3c)  and  arrive  at  the 
linearized  boundary  condition  at  the  body  surface: 


2.78 


. . 


For  first-order  axial  and  cross-flow  problems,  Eqs. 
(2.2.6c)  and  (2.2.7c)  they  take  the  form 


cnOJ.(x,R)  =  R'  (x)  (2.2.11a) 

-  <Plr(x,R)  =  -  1  (2.2.11b) 


As  remarked  by  Van  Dyke  (1951)  these  approximations  have  the 
advantage  that  the  resulting  first-order  solutions  satisfy  the 
supersonic  similarity  law  (the  supersonic  counterpart  of  the 
Goethort  rule),  which  is  not  true  if  the  exact  tangency  con¬ 
dition  is  used.  Except  possibly  for  extremely  thick  bodies, 
the  resulting  error  in  using  the  approximate  tangency  con¬ 
ditions.  Eqs .  (2 . 2 . 11)  must  be  small  at  all  Mach  numbers  (VanDyke,  lb5i). 

The  linearised  boundary  conditions  (2.2.11a)  and  (2.2.11b;  lead 
to  th'a  following  values  for  the  arbitrary  functions  f(x) 
and  gv.x;  (Von  Karman,  1932 ) i  (Tsien ,  1938): 


i  (x) 


K  ( x )  R '  f  x  j 


J_  dS(x) 
2  7T  dx 


and 

g(x)  2R(x)  K*  (x)  £  -ISOO 

7T  dX 


(2.2.12a 


(2.2. J 2b) 


2 

where  J(\)  u  R  (x)  is  the  cross-sectional  rea  of  the  body 
at  x . 
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Denoting  the  axial,  radial  and  azimuthal  perturbation 
velocity  components  by  u,v  and  w  in  this  first-order  solution, 
we  obtain 


u  =  cpQx  +  a  <plx  cos  9, 
v  =  mor  +  2  <Pircos  0* 


w  =  -  ^  sin  6 


(2.2.13a) 

(2.2.13b) 

(2.2.13c) 


where  ©Qjcp^  and  their  derivatives  are  given  by  equations 
(2.2.8)  and  (2.2.9)  with  f  (t)  and  g(t)  from  Eqs.  (2.2. 12)  . 

We  now  replace  the  straight  downstream  characteristics 
x-pr  by  y(x,r)  as  in  V/hitham's  theory  (1952).  With  this  modi¬ 
fication  u,  v  and  w  become 


u  = 


y  y 

-  J  — g-iP—  dt  +  ff-coge  J  (y-t  -f  8  r ) g '  (t)_  dt  (2.2.i4a) 
o  Y(y,t,pr)  pr  Q  Y(y,t,pr) 


y  ,  y  2  > 

1  j  (y-t  +  pr)f  (t)  dt  _  a,  cos  Q  J  (y-t+Sr)  g  (t)  dt 
r  q  Y (y , t ,pr)  Pr2  „  Y(y,t  ,r)  (2. 


Y(y,t  ,r)  (2.2.14b) 


r  _ q  sin  9  J  (y-t  +  Pr)g(t)  dtj 
pr2  Q  \  (y , t ,8r ) 


(2 . 2 . 14 .  c ) 


where 


Y(y,t,8r)  =  {  (y-t ) (y-t  +  2pr )} 


1/2 


(2.2.15c) 


y  is  now  determined  from  the  condition  that  y(x,r)  =  constant 
is  a  characteristic  curve,  that  is,  along  it 
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—  =  cot  (p.  +  F>)  , 
dr 


(2.2.16) 


where  p,  is  the  local  Mach  angle  and  6  is  the  local  flow 
\ 

direction.  The  value  of  y  on  a  characteristic  has  not  been 
defined  uniquely,  although  on  the  body  it  must  be  approxi¬ 
mately  equal  to  x  -  gr  (which  it  replaces  in  linearized  theory) 
It  is  now  made  quite  definite  by  taking  it  equal  to  the  vakie 
of  x  -  gr  at  the  point  where  the  characteristic  produced  meets 
the  body  axis. 

On  expressing  cot(u  +  5)  in  terms  of  u,v  and  w  by  means 

9  2  2 

of  Bernoulli's  equation  and  neglecting  terms  0(u"  +  v  +w  ), 
the  integration  of  eq.  (2.2.16)  with  respect  to  r  yields 


x  =  gr  -  c(y,r)  +  y 


(2.2.17) 


where 


c  (y  ,r)  -ittliS-  / 


2(3 


N(r )  -  N  (R) 

,  .,1/2 

(y-t) 


f  (t)  dt 


1/2 


,V2 


*  m2  / «  r-(r>--j^Tr2  •  N(R)+(y-t);^ 


N(r)  +  (y-t)  N  (R)  -  (y-t) 


dt 


(^+1)M4  _  f  N  (r)  -  N  (R)  ^ 

— * — ^ -  a  eos  Q  J  — - 1/2~  Svt)  dt 


2g  «  (y-t) 


j(y+l) M4  -  2M2g2 


2g 


i  rn_n  rr„h(r)  -  (y-t)1/2  N  (R)  +  (y-t ) 1/2 

V.  a  cos  0  Jini - 7  I'/? - 172 

f  lN(r)  +  (y-tr7**  N(r)  -  (y-t)A/^ 


+  —  a  cos  0  J  N (r)  (y-t)g'(t  )dt  -  ^-a  C°-S-  -  /  N(R)  (y-t  )gr(t )  dt 
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—  2M  a  cos  0 


-21 (*  CL 


cos  e 


j  ln  fsM  - 

i  »(K)-(r.)1/J, 

/-  tej ' 


^  g#-(t) 


dt 


g(t)  dt, 


(2.2. 18) 


where 


N(r)  =  N(y.,t,2pr)  =  (y-t  +  28r)1/2, 


N (R)  =  Nl  y,t, 2BE(y.)}  =  {  y-t+2BR (y )} 1/2 


(2.2.19a) 

(2.2.19b) 


and  R  =  R(y) . 


Writing 


where 


c(y,r)  =  k  J'{y,T.,R(y.),M,a,rci}  r1//2 


k  =  2  1/2  fy+i)  :m4p_3/2 


(2.2.20) 


(2.2.21) 


and  changing  the  variable  of  integration  by  substituting 


(y-t^)a/2  •,= 


(2.2.22) 


the  R- function  way  'be  written  as  -(see  .-Appendix  2.2.1) . 


K;y.,r,R(y),M,a,0}  =  -=;  [I-, (r)  -  I,{  R(y )}  ] 


[I2(r)-I2{^y)j] 


(cyrHL  )M2  ~V2Br 


-  [i2(r)  -i^RCy)}]  +  2a^y  [i2(r)  •• 


2 


F  {  y  ,r,R(y)  ,m}  .  — —  [I,  (r»  -  I,(R)j - - [I.,<r)  -  I„(R)] 

**  '  (7+1  )M2  2Pr 

(2.2.26) 


where  the  expressions  for  1^  and  Ig  are  given  by  (2.2.24a)  and 
(2.2.24b). 

If  we  assume  ftr/y  »  1,  we  obtain  Whitham’s  simplified 


F- function  which  depends  only  on  y‘. 

Vf 

F(y)  =  2  j  f'  (y-|2)d£ 


(2.2.27) 


Making  use  of  the  substitution  (22) ,  we  can  write  this  in  the 
form  given  by  Whitham  (Whitham,  1952)* ‘ 

y  . 


«y)  -  « 

„  (y-t)1/2 


(2.2.28) 


F- Function 


The  F-f unction  is  given  by  Eq.  (2.2.23)  in  which  the  I- 
integrals  of  Eqs.(2.2.24)  can  be  rewritten  after  substitution 
of  the  values  of  f(x)  and  g(x)  from  Eqs.  (2.2.12) 


In  (r)  =  —  J  (42-K2pr)1/2  s"(y-£2)  d£, 

^  2  ir  , 


(2.2.29a) 


I2(r)=  -i-  /  (|2  +  23r)  VV  (y-£2)  d£ , 

^  277 


(2.2.29b) 


J 

I1(r)=  -2  /  (f2  +  23r)1/2s"(y-^2)  d£  =  2(3  I^r),  (2.2.29c) 


I„ (r)  =  -  /  (£2 + 2pr)"1/2s'  (y-£2).  d£  =  20I2(r),  (2.2.29d) 
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.1/2 


I3(r) 


a 

7T 


/  V 


i^r)1/2 


(y-C2) 


where  primes  indicate  differentiation  with  respect  to  the 
2 

argument  (y-£  ). 

The  I-(R),  Ip(R)  etc.  are  obtained  from  above  by  replacing 
r  by  R(y) . 

If  f(x)  and  g(x),  which  are  proportional  to  S^x),  are  given 
as  polynomials  we  proceed  with  the  integration  as  given  below. 
If,  however,  the  value  of  (x) ,  which  is  assumed  to  be  con¬ 
tinuous,  is  known  for  at  least  (n+1)  values  of  x,  say 
xo,xl’  xn’  the  simPlest  and  most  often  used  process 

consists  in  selecting  an  approximate  function  yR(x)  which 
takes  on  the  same  values  as  does  S* (x)  for  each  of  those 
(n+1)  values  of  x.  Here  again  the  choice  of  polynomials, 
called  the  collocation  polynomials,  is  convenient.  Fox, 
whereas  in  the  general  case  there  may  be  no  function,  or 
there  may  be  several,  the  existence  and  uniqueness  theorem 
states  that  there  is  one  and  only  one  polynomica 1  of  degree 
n  ox-  less  which  takes  on  the  proscribed  values  at  each  of  the 
(n+1)  points.  The  basic  interpolation  formula  was  given  by 
Lagrange  which  expresses  the  approximate  polynomial  in  terms 
of  only  the  ordinates  given.  It  is,  however,  often  con¬ 
venient  to  use  Newton'S  interpolation  polynomial  with 
divided  differences  or  Aitken'S  iterated  interpolation 
polynomials.  These  latter  formulas  are  especially  suitable 
for  the  detection  and  propagation  of  error  at  each  stage. 
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For  details  reference  (Hildebrand , 1956)  may  be  consulted.  When 
the  values  as  well  as  the  derivatives  upto  a  certain  order  of 
a  function  are  given  at  specified  arguments,  we  can  use 

Hermite’s  interpolation  polynomials  (also  called  osculating 

✓ 

pblynomials)  to  approximate  the  given  meridian  curve. 

We,  therefore,  write 


s'(y-£2)  =  s'{(y  +  2£r)-  (£2  +  2£r£ 


N 


-z 


An(y,£r)  (£2  +  23r)n 

(2.2.30) 


Substituting  this  or  its  derivative  in  the  expressions  (2.2-29) 

it  is  readily  seen  that  all  the  integrals  of  I-^Cr),  I2(r), 

I^(r)  and  I2(r)  reduce  to  the  following  form; 

1/2 

y 

Tp/g (t)  =  /  (£2  +  2£r)p/2d£,  p=-  1,0, 1,2,  .  .  .  ,N  (2.2.31) 
o 


The  reccurence  formula  for  the  ~  integrals  is  easily 

established  as 


T  /2(r>  -  Z 

P/2  D  +  1 


26rp 


(2.2.32) 


and  since 


1/2 

y 

T_l/2  (r)  =  /  (e2  +  2Pr)"1/2de  =  tn 


+  (y  -f  2gr) 1/2 
(2£r)1/2 


(2.2.32) 


all  the  integrals  of  I1(r),  I2(r),  I^(r)  and  I2(r)  can  be 

written  in  closed  form. 

The  integrals  of  I2(r)  will  be  of  the  form 
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!  e3(?2+  26r)p/2  di 


(2.2.34) 


which  is  readily  integrable  as  it  can  be  reduced  to  a  standard 


form. 


Ho.  i  2  +  2 

i  /2  =  (y  +  2{3r)2  ‘  %  - 

P/^  (p-f-2)  (p+4)  (p+2)  (p*4) 


(2.2.35) 


The  F- function  of  Eq. (2.2.23)  can,  therefore  be  determined  in 
all  cases  where  the  body  is  smooth,  in  closed  form. 

The  F-curve  gives  immediately  a  rough  description  of 
the  flow  pattern  since  it  shews  whether  the  characteristics 
are  converging  in  compression  (dF/Sy  >  0)  then-  a  shock  will 
appear  or  diverging  in  expansion  (dF/dy  <  0). .  The  method  of 
area  balancing  given  by  Whitham  (1952)  can  be  used  to  build  up 
the  details  of  the  flow  pattern. 


The  Pressure  Signature 


From  Bernoulli’s  equation,  which  is  approximately  true 
since  the  small  entropy  changes  at  a  shock  contribute  a 


term  of  smaller  order, 


_X_ 
2  -,-l 


42  =  r/l  +  1^(1  -  ^5)}  1  - 

n  w  n  tt* 


(2.2.36) 


where  Ap  is  the  pressure  in  excess  of  the  undisturbed  pressure. 
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(1-2  a£  cos  B) 


A£  =  (2&)1/4y 
P0  (l+7> 


r-3/4 

1/2 


/  F(y)  dy 


1/2 


(2.2.39) 


where  yQ  is  the  first  zero  of  F(y) ,  the  Whitham  function  for 
far-field  points.  In  between  the  shocks  equation  (2.2.38)  holds. 
The  two  integrals  of  Eq.  (2.2.38)  may  be  evaluated  in  closed  form 
using  the  T- integrals  of  Eqs.  (2.2.32)  and  (2.2.33). 


Examples 

In  order  to  compare  the  near-field  F-curve  with  the  far- 
field  Whitham  function,  we  have  calculated  these  functions  and 
the  pressure  signatures  for  two  examples  of  body  of  revolution, 
(a)  The  shape  of  the  first  body  is  given  by 

2 

R(x)  =0.1  (x-  ::  ) ,  0  <£  x  ^  1  (2.2.40) 

=  0  ,  x  >  1 

The  Whitham  F-function  is  easily  calculated  analytically  and 
is  given  by 


F(y)  =  0.02  [  (1  -  6y +  6y2)  yy2  + 2(l-2y)y3/2+ 1.2  y5/2] 

(2.2.41) 

The  F-function,  as  given  by  Eq.  (2.2.23)  has  been  expressed  in 
terms  of  the  1^,  Ig  and  1^  integrals  which  are  in  turn  ex¬ 
pressible  by  means  of  the  T  and  Q  integrals  of  Eqs.  (2.2.32) 
and  (2.2.34).  For  this  example,  if  we  put 

ar  =  y  +  2£r  and  aj,  =  y  +  26R(y)  (2.2.42) 
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=  l  +  2(m  +a  cos0  +  a$  cos  0 )  +  higher  order  terras 

ox  XX  or 

and 


1  -  -  -  2  O  +  at<p  x  cos  0  +  a  <por  cose] 

u 

neglecting  higher  order  terms. 

Hence,  within  the  order  of  accuracy  employed  in  the 
first-order  theory  equation  (2.2.36)  may  be  approximated  by 

^  =  -  7M2  [{0ox  +  (<p^x  +  ®or^aCOS®]  (2.2.37) 

Po 

Substituting  the  values  of  the  derivatives  of  cd0  and  cpj  from 
Eq'2.2.9)  and  making  use  of  (2.2.12)  and  (2.2.22),Eq.  (2.2  37)  may  be 
rewritten  as 


I 


1/2 

s" (y-£2) 
(#-2+2pr) 


1/2 


d*  - 


3q  cos  9 

,5r 


1/2 

/  „e| 

„  (f2+2pr)1/2 


(2.2.38) 


Since  the  points  of  a  iront  shock  are  "effectively  at  large 
distances",  Whitham’s  approximate  treatment  (1952)  is  valid  on 
the  shock  and  we  reproduce  his  result  with  minor  alteration 
for  the  inclined  body.  Hence  for  the  front  shock,  when  r  is 
large 
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For  the  near  field,  if  we  put 


A  =  1  -  a  ,  A_  =  1  -  aD  (2.2.47) 

r  r  '  R  R 

where  a_  and  a^  are  given  by  eq.  (2. 2.42)  the  analytical  expression 
for  the  integrals  required  for  the  calculation  of  the 
F-curve  are  as  follows. 

I2(r)  =  0.01  [Ar(15A^- 6)T1/2(r)  +  (60A^-6)T3/2(r) +  90A^T5/2(r) 

+  60ArT?/2(r)  -f  15Tg/2(r)J,  (2.2.47a) 

I2(r)  =0.01  [3A^(1- A^)T_1/2(r) +3Ar(2- 5A^)T1/2(r) +3(l-10A^)T3/2 

-  30  A^/2<r)  "  15ArT7/2(r)  -  3T9/2(r)] 

(2.2.47b) 

I3(r)  =0.01  [Ar(15A^-6)Qv2(r)  +  (60A^-6)Q3/2(r) +90A^Q5/2(r) 

+  60ArQ?/2(r)  +  15Q9/2(r)l  (2.2.47c) 

Agkin  the  values  of  I1{R(y)},  I2{R(y)}  and  I3{  R(y)}  are  ob¬ 
tained  from  (2.2.47)  by  replacing  r  by  R(y). 

The  pressure  signature  in  between  the  shocks  can  be 
expressed  in  terms  of  the  known  integrals  given  above 
yielding 


Ix(r)  =0.01  [(l-6ar  +  6aJ)T1/2(r)  +  (6-12ar)T3/2(r) +6T5/J;i0j, 

(2.2.43a) 

I2(r)  =0.01  [(ar-3a2  +  2a^)T_1/2(r)  -  (l-6ar  +  6a2)T1/2(r) 

(3  -  6ar )Tg^2  (r )  -  2  T5/,2(r)], 

(2.2.43b) 

I3(r)  =  0.01  [(l-6ar  +  6a2)Qi/2(r) +6(l-ar)Q3/2(r) +6Q5/2(r)] 

(2.2.43c) 

The  I-j{R(y)},  I2{R(y)}  and  Ig{R(y)}  are  obtained  from  (2.2.43) 
by  replacing  r  by  R(y). 

The  pressure  signature  in  between  the  shocks  is  given  by 


^  =  0.02  7M2[{1+  (3p)3COS0}{  (l-6ar  +  6a^)  T_1  2(r) 

P0  •  * 

+  (6-12ar)T]/2(r) +6T3/2(r)}  -  |  a  cos  0  Ij(r)J  (2.2.44) 


For  the  front  shock  the  jump  in  pressure  can  be  calculated 
from  eq.  (2.2.39). 

(b)  The  shape  of  the  second  body  considered  is  given 
by 


R(x)  =  0.1  [1  -  (1-x)  ],  0  ^  x  ^  1 


(2.2.45) 


=  0.1 


x  >  1 


The  Whitham  F-function  for  the  far-field  points  is  given  by 


F(y)  =  [0.3(l-y)4-0.l2(l-y)]y1/2  +  [0.4 (1-y)3- 0. 04j  y3/2+ 

+  0.36 (l-y)2y 5/2  +0.171(l-y)y7/2+  ~  y9/2  (2.2.46) 

Ou 
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The  pressure  signature  in  between  the  shocks  is  given  by 

££  =  0.027M2  [{  1  +  (33)  a  cos  e}{  (15a  J  -  30a2  +  90a  J  -  54ap  +  9)T_1/2 
po 

-  (60a2  -  180a2  +  lSOa^  +  54)  T1/2(r) 

+  (SOa^  -  180ar  +  90)T3/2 (r ) +  60 (l-ar )T&/2 (r)} 

-  ~  a  cos  e  I1(r)]  (2.2.48) 


Conclusions 


The  worked  examples,  given  in  the  paper,  show  a  signi¬ 
ficant  difference  between  F-function  in  near  and  far  field 
points.  It  is  quite  clear  from  the  plots  that  for  0  =  180° 
at  a  certain  angle  of  attack,  the  peak  pressures  for  near 
fields  are  higher  than  those  at  zero  angle  of  attack.  From 
the  plot  of  F-function  and  also  from  the  pressure  plots,  it 
is  seen  that  the  expansion  is  stronger  than  the  case  of  zero 
angle  of  attack.  This  is  even  otherwise  obvious  from  the 
physical  point  of  view.  For  9  =0°,  the  reverse  is  true, 
i.e.,  the  pressure  peaks  become  smaller  and  expansion  becomes 
\  iker  with  angle  of  attack.  Thus  we  observe  that  with  an 
angle  of  attack  we  have  a  shock-expansion  system,  in  which 
the  strength  varies  as  the  azimuthal  angle  is  varied  around 
the  body  of  revolution. 


(r) 
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2.2.4;  Pressure  signature  for  the  body  R(y) 
with  and  without  angle  of  attack  (r 
all  cases) . 


figure  2.2.6:  Pressure  signature  for  the  body 
R(y)  =  0. 1[  l-(l-y)J]  with  and  without,  angle  of 
attack  (r  =  0.125  for  all  cases). 


2.2.2  DETERMINATION  OF  A  BODY  OF  REVOLUTION  FOR  A  SPECIFIED 
PRESSURE  SIGNATURE  AT  ANY  FIELD 

Summary 

Design  of  bodies  of  revolution  to  produce  specified  sonic- 
boom  signatures  at  a  given  Mach  number  and  lateral  spacing  from 
the  body  uses  the  far-field  theory  of  Whitham  and  a  procedure 
has  been  put  forward  by  Barger  (1968)  for  calculating  the 
shape  of  the  bodies  of  revolution.  Since  the  above  theory  does 
not  strictly  apply  to  near-field  region,  the  present  work  gives 
a  mathematical  procedure  for  determining  the  bodies  of  revolution 
for  specified  S' ^natures  at  any  field.  However,  it  should  be 
pointed  out  here  that  the  near-  and  far-field  theoreis  merge  at 
a  certain  spacing  from  the  body  and  beyond  this  field  the  far- 
field  theory  due  to  Barger  could  be  used  for  its  mathematical 
simplification. 
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Introduction 

it 

The  Whit  ham  Theory  (Whit  ham,  1952)  has  been  shown  t  •>re-.i 
diet  accurately  far-field  signatures  produced  by  a  body  of  ^ 
revolution.  This  theory,  however,  is  not  very  accurate  in  the 
near  and  mid-field  points  and  Sec.  2.2.1  points  out  the  limitations 
and  applicability  of  Whitham' s  far-field  theory.  The  far-field 
theory  applicable  to  the  inverse  problem  for  determining  the  . 
shape  of  a  body  of  revolution  that  will  produce  a  specified  .• 
pressure  signature  at  a  given  Mach  number  and  certain  field 
away  from  the  body,  has  been  given  by  Barger  (1968).  This  theory 

l  .  ** 

does  not  strictly  apply  to  near-  and  mid-field  points.  The 

present  work  employs  near-  and  mid-field  theory  in  order  to 

>  '  . 

determine  a  body  of  revolution  for  a  specified. physically  . 

j  n 

obtainable  pressure  signature  at  the  near-field  region.  Although 
variations  in  the  shape  of  a  generating  body  do  not  normally  : 
prevent  its  signature  from  developing  into  the  well-known  N-shaped 
wave  at  large  distances,  there  are  some  signatures  that  do  not 
attain  the  classical  N-wave  form  in  the  far-field.  One  such 
example  is  given  in  the  reference  (Barger,  1968)  where  the  specified 
near-field  pressure  signature  has  a  form  with  its  positive  and 
negative  parts  separated  by  a  section  of  zero  overpressure. 

Moreover,  even  if  the  signature  does  approach  an  N-wave,  it  may 
develop  so  gradually  that  it  does  not  attain  its  final  asymptotic 
form  within  the  distance  corresponding  to  objectionable  overpressures 
(McLean,  1965),  or  even  if  N-wave  io  formed  at  the  given  distance 
the  overpressures  are  of  less  magnitude  than  annoying  levels. 
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Tliereis  a  need  for  analytical  methods  that  are  capable  of 
accurately  predicting  the  near  flow  field  pressures  about 
arbitrary  aircraft  configurations  for  design  optimizations  in 
connection  with  sonic  boom  as  well  as  the  performance  of  the 
vehicle  (Hunton,  1968).  Near-field  pressure  signatures  are  also 
of  interest  in  connection  with  certain  sonic-boom  studies  such 
as  wind-tunnel  investigations  and  signature  measurements  using  a 
probe  aircraft.  Available  theories  are  either  too  cumbersome 
mathematically  for  ease  of  handling  on  the  computer  (e.g., 
method  of  characteristics)  or  are  limited  in  range  of  appli¬ 
cation  (e.g. ,  Whitham,  1952)  because  of  some  mathematical  simpli¬ 
fication  such  as  the  linear  theory  concept.  As  a  result,  problems 
of  analysis  do  occur  in  the  flow  regions  very  near  the  aircraft 
in  connection  with  the  design  and  arrangement  of  configuration 
components.  Various  wing- body  combinations  have  to  be  examined 
to  come  up  with  various  desired  configurations.  Whitham’ s  theory 
is  easily  applied  to  the  bodies  of  revolution  but  can  not  be 
applied  to  wings  to  determine  the  flow  fields  at  the  near-field 
points.  For  these  fields  a  wing  can  not  be  replaced  by  an  equi¬ 
valent  body  of  revolution  and  thus  new  prediction  techniques 
wherein  a  wing  is  represented  by  some  sort  of  spatial  distribution 
of  singularities  over  the  entire  wing  surface,  have  to  be  developed 
This  phase  of  the  sonic  boom  research  is  now  under  progress. 

Thus  the  present  work  is  a  preliminary  study  towards  the  overall 
optimization  of  the  aix’craft  in  connection  with  sonic  boom. 

The  purpose  of  the  present  paper  is  to  describe  a  procedure 
for  designing  bodies  of  revolution  corresponding  to  specified 
pressure  signatures  at  any  field.  "* 


Let  the  steady  stream  have  velocity-  U  in  the  x-direction. 


and  at  a  general  point-(x,r)  let  the  velocity  be  (U+Uu,Uv). 
Assuming  the  flow  to  be  irrotational  the  perturbation  velocities 
u  and  v  may  be  deduced  from  a  potential  0  which,  on  the  linear¬ 
ized  axisymmetric  theory  satisfies  the  equation 


+  ;  rfr  '  '  0  (2.2.19, 

where  3^  !  =  (M2  -  1)  and  suffixes  denote  partial  differentiation. 
The  solution  of  Eq. (2.2.49)  which  represents  a  disturbance 
propagated  downstream  from  a  body  is 


x-0r 


f(t)  dt 


<x-t)2  -  p2r2}1/2 


(2.2.50) 


giving 


x-j3r 


u 


f  (t)  dt 


{(x-t)2'  -  r-2r2}1/2’ 


(2.2.51) 


x-csr 


v  = 


-I 

r 


(x-t)  f  (t)  dt 

S  (x-t)2  -  A2; 


(2.2.52) 


The  downstream  approximate  characteristics  of  the  equation  ax*e 
x  -  0r  =  constant .  We  now  replace  x-|3r  by  y  (x,r),  a  better  * 
approximation  to  the  exact  characteristic.  Carrying  out  the 
modification,  u  and  y  become 


u 


--  / 


f  (t)  dt 


[  (y-t)  (y-t  +2pr>] 


TTZ- 


(2.2.53) 


y 

1  (y-t  +  pr)f  (t)  dt 
r  Q  [  (y-t)  (y-t  +  2fir) ) 1/2 


(2.2.54) 


y  is  now  determined  from  the  condition  that  y(x,r)  =  constant 
is  a  characteristic  cusxo,  that  is,  along  it  dx/dr  =  cot  (p.+e), 
where "p. -is  the  local  Mach  angle  and  9  is  the  local  direction  of 
flow.  Proceeding  as  detailed  in  Whit ham* s  paper  (Whit ham,  1952) 
and  assuming  pr/y  »  1,  u,v  and  the  equation  of  the  characteristic 
are  given  by 


U  ss 


.  -  JM 


(2p) 


1/2 


pu. 


(2.2.55) 


l 

x  =  0r  -  k  F(y)  ra  +  y. 


(2.2.56) 


where 


I 

y, 

F(y)  ='  /  (2.2.57) 


and 


k  «  2~1/2(tKL)mV'3/2 


(2.2.58) 


From  Bernoulli's  equation,  which  is  approximately  true  since 


the  small  entropy  changes  at  a  shock  contribute  a  term  of  smaller 

2  2 

order,  we  have  by  neglecting  terms  0(u  +  v  ),  2 


AH  =  ,-y 


(2.2.59) 


where  Ap  =  p  -  pQ  is  the  pressure  in  excess  of  the  undisturbed  ' 
pressure  pQ. 

Then  the  problem  is  to  solve  the  singular  integral  equation 

s  ’  ,  -  r 

(2.2.53)  when  u(y)  is  given  by  (2.2.59). 

To  Find  a  Body  of  Revolution  for  a  Given  Near-  or  Hid- Fie Id  Pressure 
Signature’* 


Solution  of  the  Integral  Equation 


Method  1  r~  Analytical*  -  Numerical  Method: 


The  integral  equation  now  to  be  solved  is 


_llfvv  _  f  f  ft)  dt _ _ 

uvy;  —  J  •  *  ---  1/2  l/lt 

n  (y-t  1 2pr  )L/Z  (y-t ) 


(2.2.60) 


where  u (y)  is  known  from  Eq.  (2.2.59)  and  £  (y)  is  the  unknown  function. 
Putting  T(t)  =  -£*  (t)  for  convenience,  Eq.  (2.2.60)  may  be 
rewritten  as  - 


u(y).  =  /  K(y.,t)  f(t)  dt. 


(2.2.61) 


in  which  the.  kernel  has  the  form 


. ««*' . . . * . . . i . . . . . I . . . . . . . . . . . . I" 


(2.2.62) 


K(y,t)  =  -.p-(y-j-/2, 
(y-t)  ' 


where 


G  (y ,  t )  =  (y-t  +  2pr)“iy2 


(2.2.63) 


is  continuous  in  T(0  ^  t  ^  y  ^  b(say)). 

Precisely  as  in  the  case  of  Abel’s  equation,  of  which  Eq. 
(2.2.61)  is  an  immediate  generalization,  we  see  that  Eq. (2.2.61) 
cannot  have  a  continuous  solution  unless  u(y)  is  continuous  in 
1(0  £  y  £  b)  and  u(0)  =0.  We  will  assume  that  these  conditions 


are  satisfied  by  u(y). 


Multiply  Eq. (2.2.61)  by  (z-y)iy2,  where 


0  <  z  ^  ba 


and  integrate  the  resulting  equation  with  respect  to  y  from  0 


to  z,  obtaining 

z  z 

r  U  (y )  dy  _  f  1 

;  cz-y)1/2  =  { 


■  ^V/2  dt^y- 

(y-t)1/2 


(2.2.64) 


The  right-hand  side  of  this  formula  reduces  by  Dirichlet’s 
extended  formula  (Tricomi,  1957  )  to 


/  *«>  / - V/2’t)  l/z  “ydt 

„  t  (i=-y)  /  (y-t)  / 


If  we  then  write 


(2.2.65a) 


F(z)  =  J  --atffa,  dy, 

o  <z-y> V 


z 

L(x,t)  =  / - - T/7  dy  (z>t) 

t  (z-y)1/Z(y-t)1/Z 


equation  (2.2.64)  takes  the  form 


F(zO  =  /  L(z,t)  f(t)  dt 


(2.2.65b) 


(2.2.66) 


We  will  now  show#  that  the  kernel  L  of:  this  equation  is  con¬ 
tinuous  in  T,  except  possibly  on  the  line:  z.  =  t  where,  it  is 
not  yet  defined.  For  this  purpose  introduce  £  as  the  variable 
of  integration  in  place  of  t  by  means  of  the  formula 


z-t 


(2.2.67) 


We  thus  get  when  z  >  t 


L(z,t)  =  J 


r  G[.  (z-t)  4+t,t 


Ur  d  ^ 


(2^2.68) 


Since,  this  integral  remains  convergent  when  we.  replace  G  by  the 
upper- limit,  of  its  absolute,  value,  it  follows  that  it  is  uni¬ 
formly.  convergent  in  T,  and  since  G  is  continuous,  L  is  also 
continuous  wherever  it  is  defined;- 

Iii  order  next  to  see  whether  L  approaches:  a.  limit  as  the 
point;  (z,t)  approaches  a  point;  (c,c)  on  the  hypotenuse  of  T, 
we  apply  the  mean  value  theorem: for  integrals  to  the  original 
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definition  of  L  given  in  Eq. (2.2.65b) ,  getting 


L(z,t)  =  G  (x,t ) 


z 

{  (z-y ) (y-t ) 1/2  “  sin  v/2 


G(s,t) ,  (t<s<z). 

(2.2.69) 


Consequently 


lim  L  (z ,  t )  =  it  G  (c ,  c  )  =  — - 

z->c  (23r)v 

t->c 


(2.2.70) 


from  equation  (2.2.63). 

We  see  then  that  L  is  continuous  throughout  T,  and  that 
if  we  have  G(y,y)/0  at  any  point  of  I,  it  follows  that 
L(y,y)  ji  0  in  I. 

In  order  that  the  equation  (2.2.63)  be  a  Volterra  integral 
equation  of  the  first  kind  with  finite  kernel,  G  should  be  such 
that  L(z,t)  have  a  partial  derivative  with  respect  to  z  finite 


in  DCo^t^b^  and  whose  discontinuities, 
distirbuted. 


if  any,  are  regularly 


Let 


G1  (z,t)  =  —  =  -  -  (z-t*2f3r)-3/2 
1  dz  2 


(2.2.71) 


which  is  obviously  continuous  in  T.  If  we  now  differentiate 
Eq.  (2.2.68)  with  respect  to  z  tinder  the  integral  sign,  we  get 


/  (— ^t) 

o 


1/2 

G1[(z-t)  4+t,tj  d£ . 


(2.2.72) 
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Since  this  integral  is  obviously  uniformly  convergent,  we  see, 
by  the  ordinary  test  for  differentiating  an  infinite  integral 
under  the  sign  of  integration,  that 


Lj  (z ,  t ) 


r  t  V2 

—  =  J  <-*-)  G.ffz-t)  £+t,t]  d£, 

o 


or 


JL^(z,t)  —  -  —  J  ) 


1/2 


_dl 


2  1-e  {  (z-t)  £+2£r}3/2 


which  is  valid  throughout  T. 

It  can  now  be  shown  that  the  continuous  solution  of  Eq. 
(2.2.66)  satisfies  Eq.  (2.2.61)  and  that  it  is  unique  (Bochei; 

1909) . 

We  shall  now  reduce  equation  (2.2.61)  to  Volterra’s- integral 

I  ’ 

equation  of  the  second  hind  which  is  more  readily  solved  by 

* 

Picard’s  method  or  numerically. 

Differentiating  Eq. (2.2.66)  with  respect  to  z,  we  have 

z 

F'(z)  =  L(z,z)  ¥(z)  +  J"  (z , t )  ¥(t)  dt  (2.2.74) 

o 

But  by  Eq.  (2.2.70) 


i»  (z ,  z )  =  1 75  /  0 

(.2$r)1/d 


(2.2,75a) 


It  can  be  easily  shown  that  (Badaaard,  1928) 
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M 

°  (y-t)* 


dt  „  jmj-zMa  dt  +  F(y)  /— dt 


(y-t)* 


-+^/s 


(y-t) 


Using  the  above  formula  we  get, 


Pfe)  =  f  /  dy  =  -  i  /aizL0 

dz  o  (z-y)1'  2  o  (z-y)d/ 


1  f  u(y)  -  u  (z  ) 

2  i  (z-y)3/2 


Hence  Eq. (2.2.74)  becomes 


dy  + 


ufo) 


(2.2.75b) 


P  (z)  = 


(20r)' 


(z)  +  /  L-,  ( z,t)  ¥(t)  dt 


(2.2.76) 


where  F  (z) ,  the  known  function  is  given  by  Eq. (2.2.75b)  and 
the  kernel  L^s.t)  by  Eq.  (2.2.73) 

The  Kernel  Function -L^ (zst) 


L.  (z,t)  =  -  -  / 


2  o  (l-£)1/2{  (z-t)  £+20r} 3/2 


(2.2.77) 


Substituting  £  =  sin  e  and  simplifying,  we  have  for  z>t  (see 
Appendix  2.2.2) 


Lj (z , t )  = 


(z-t)  (z-t  +  2pr) 


T75  E(k) - r/2  K(#c) 

1/2  (z-t)  (z-t  +  2£r)1/2 

(2.2.78) 


where 
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I  W 


K  =  ( - )  (2.2.79/ 

z-t+2pr 


and  K(/c)  and  E(k)  are  complete  elliptic  integral  of  the  first  , 
•*nd  the  second  *kind. 

It  remains  to  find  the  value  of  (z , t )  for  z=t.  We  have 

t 

from  Eq .  (2 . 2 . 73  ) 

4(2  pr) 3/2 

(2..2.80) 


Method  2  -  Direct  Numerical  Method  of  Solution 

(Collocation  Method) 


;[L,(z,t)]  =  -  ~  - - 

1  z-t  2  (2pr)  - 


J  <A-¥2  dt  =  - 

■o  i-4 


It  is,  however,  possible  to  solve  the  integral  equation  iby 
application  of  the  collocation  method  ‘Without  reduction  to  the 
second  kind. 

The  integral  equation  to  be  solved  is 


S 


o 


f'  (t  )  dt  _ 

Cy-t+2(3r)1/2(y-t)1/2 


=  -  u(y) 


(2.2.81) 


where  u(y)  is  assumed  to  be  a  'known  continuous  function  in 
ItO^y^b(say))  and  u(o)  =  0. 

In  the  general  case  we  can  assume  that  f '  (t)  can  be  ex¬ 
pressed  in  the  form 


f'(t)  =  j 


(2.2.82) 


which'  includes  polynomials  as  special  cases.  We  proceed  to 
determine  the  unknown  coefficients ,  an  by  introducing  expression 
Eq. (2.2.82)  into  Eq. (2.2.53),  giving 


«  y  2 

/  an  J  T/2  T/?  52  "  b  (y )  • 

£o  o  (y-t+20r)1/2(y-t)1/- 


(2.2.83) 


The  following  two  types  of  integral!  occur  in  Eq. (2. 2, ^3) 
according  as  n  =  2m,  an  even  positive  integer  or  n  =  2m+I.,  an 
odd  positive  integer, 


I,  (y)  =  / 


tmdt 


(y-t+2pr ) 1/2 (y-t ) 1/2 ’ 


(2.2.84a) 


2m+l 


Io(y)  =  J  - 1 - TTo - 772  *  (2.2.84b) 

o  (y-t+2pr  ) 1/2  (y-t ) 1/2 

In  Appendix  2.2.3,  it  is  shown  that 

I1(y)  =  2[y“To-  (“)  ...  +  <-l)r(*)  Y®“rTr  ...  +(-l)®Tmj 

(2.2.85) 


where 


1/2  1/2 
in  (^ - =75 — ), 


and  the  recurrence  formula  for  the  evaluation  of  T,  is 

M 

2m- 1 

(ZgKU. 


(2.2.86a) 


<2. 2. 86b) 


Y  =  y  +  2£r 


(2; 2. 86c) 


and 


2  m+1 

l2<y>  =  <y+2Pr>1/2  A2b+2 


where 


2bH-2 


2n(l+k1)  A2g- 

(2m+l)  k? 


with 


A2  =  \  [K(kx)  -  ECkj)], 


A4  =  -ij  f  (2+k2)  K(kj)  -  2 (1+k2)  ECkj)] 
3k  j 


in  which  the  modulus  k1  is  given  by 


(2.2.87) 


(2.2.88a) 


(2.2.88b) 


(2.2.88c) 


k?  =  — ^ —  ,  0  <  <  1  (2.2.88d) 

y+2pr 

and  K(k1)  and  Efkj)  are  complete  elliptic  integrals  of  the  first 
and  second  kind. 

The  integral  equation  (2.2.53),  then  reduces  to 
N 

5  an^n(y)  =  "  u(y) 


(2.2.89) 


where  Jin\y)  arc  Known  for  all  values  of  y  in  I.  Chosing  N 
suitable  values  of  y,  (y-^ygj  yN>  for  which  u(y)  is  known, 

we  have  the  system  of  equations  to  determine  an, 


y  a2'/^'2^y2^  + 


.  +  a. 


+  aNu 


4(y2> 


-  uCyp, 


=  "  u  (y0) , 


af^lCy^)  +  a2^2^N^  + 


■  al!4(yH) 


-  u(yN), 


which  can  be  written  in  matrix  notation  as 


[Jh  {a  )  *  ~  {u  } 


{a}  -  -  } 


(2. 2. SO) 


where 


[A]  =  [„iCys)]:  is  a  (NxN)  square  matrix  of  known  elements 

{u}  =  .fu(y  )}  :  is  a  (Nxl)  column  matrix  of  known  elements 
s 


{ai  =  {ar} 


is  a  (Nxl)  column  matrix  of  the  unknown 
coefficients. 


Hence  f  (y)  is  approximately  determined,  from  which  the 
geometry  of  the  equivalent  body  of  revolution  may  be  found. 
The  use  of  the  above  method  has  been  illustrated  in  the 
"Numerical  Application"  later  in  the  section. 
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Numerical  Application 


la  order  to  illustrate  the  methods  of  solution  for  obtain- 
ing  the  shape  of  a  body  of  revolution  that  will  generate  a 
specified  signature  at  a  given  Mach  number  and  lateral  spacing 
from  the  body  the  following  example  is  chosen. 

Two  of  the  specified  pressure  signatures  as  given  by.  Barger 
(196,8)  have  been  selected  -  one  with  a  finite  rise  and  thp 
other  with  a  plateau  with  its  peak  same  as  that  for  the  finite 
rise.  The  sole  reason  for  selecting  these  pressure  signatures 
at  M  =  1.2  and  r/t  =4.2  is  that  with  these  conditions,  the 
specified  and  experimentally  measured  pressure  signatures  agree 
very  well  with  some  rounding  off  of  the  corners. 

These  signatures  have  a  stretch  of  15.3  cms  starting  from 
zero  strength  at  x-  fir  =  0  and  coming  back  to  zero  strength  at 
x  -  0r  =  15 . 3  cms .  The  peaks  of  the  signatures  are  obtained  by 
joining  straight  lines  rather  than  rounding  them  at  the  peaks 
(Barger,  1968).  This  is  done  only  for  mathematical  convenience. 
The  desired  finite  rise  pressure  signature  at  r/t  =  4.2  and  M  =  1.2 
is  given  by, 


AP  _  0.0096 


0  ^  x  £  4  cms. 


=  0.0096  -  --0-192  (x-4 ) ,  4  ^  x  £  11.5  cms. 


=  -  0.0096  +  (x-  11.5),  11.5^x  15.3  cms. 


and  the  desired  plateau  pressure  signature  at  r/l  =  4.2,  M  =  1.2 
is  given  by, 


£p  _  0 . 0096 
PQ  "  0.5 

=  0.0096 

=  0.0096 

=  -  0.0096 


x  ,  0.5  ems 


0.5  ^  x  £  4  ems 


0. 0152 
7.5 

0,0096 

3.8 


(x-4),  4  £  x  <£  11.5  ems 

(x~11.5),  11.5  £  x  £  15.3  ems 


These  two  signatures  are  plotted  in  Fig.  2.2.7.  Since  at  these 
conditions  the  far-field  theory  given  by  Barger  (1968)  is  adequate 
for  determining  the  shape  of  the  body  of  revolution,  it  is  of 
no  use  to  modify  this  far-field  theory  at  this  field  r/l  =4.2 
with  M  =  1.2  or  any  field  away  from  r/l  -  4.2.  The  body  obtained 
by  use  of  the  far-field  theory  does  not  give  any  significant 
difference  between  the  specified  pressure  signature  and  pressure 
signature  obtained  by  use  of  the  near-field  theory  (Sec.  2.2.1), 
or  in  other  words,  the  far-field  and  near-field  merge  at  r/l  =4.2 
for  this  body  at  the  given  conditions.  So  the  near-field  theory 
is  of  use  only  for  fields  between  r/l  =  0  and  r/l  =  4.2  at 
M  =  1.2.  The  numerical  procedure  adopted  to  illustrate  this  is 
as  follows: 

Barger's  procedure  is  used  to  determine  the  shape  of  the 
body  of  revolution  by  the  use  of  his  far-field  theory.  This 


procedure  involves  the  construction  of  F-f unction  and  then  the 
F-function  integral  is  inverted  to  give  the  shape  of  the  body 
of  revolution.  Then  the  Whitham's  far-field  theory  (Whitham, 

1952)  is  used  to  compute  the  pressure  signature  of  this  body  at 
r/i  =  0.5  and  also  the  near-field  theory  (Sec.  2.2.1)  is  used 
to  compute  the  same  at  r /l  =  0.5.  The  two  pressure  signatures 
are  plotted  in  Figure  2.2.8*  It  is  seen  that  there  is  signi¬ 
ficant  difference  between  the  two  signatures.  This  suggests 
that  the  far-field  theory  is  not  accurate  enough  in  these  near¬ 
field  points.  Of  course,  here  the  body  is  approx:  .lately  8-10% 
thick  and  thus  it  is  expected  that  thicker  bodies  (about  20% 
thick)  will  give  larger  errors  between  the  pressure  signatures 
obtained  by  the  use  of iar-field  and  near-field  theories. 

Now  the  far-field  pressure  signature  at  r/l  =  0.5,  M  =  1.2 
becomes  the  specified  pressure  signature  in  that  field.  We 
can  use  "method  of  collocation"  or  "direct  numerical  method"  to 
determine  the  shape  of  the  body  of  revolution.  Since  in  the 
near-field  points  the  F-function  and  pressure  signature  have  no 
simple  relationship  as  in  the  case  of  far--field  points,  the 
integral  equation  (2.2.60)  is  solved  assuming  that  the  character¬ 
istics  are  straight.  The  body  obtained  from  this  procedure  is 
the  first  approximation  and  then  the  near-field  theory  (Sec.  2.2.1) 
is  used  to  calculate  the  pressure  signature  at  the  above  field 
r/t  =  0.5  at  M  =  1.2.  It  is  found  that  this  signature  and  the 
specified  signature  do  not  coincide,  thus  suggesting  that  the 
characteristics  are  not  straight  at  this  field.  Although  the 
peaks  of  the  signatures  are  the  same,  there  is  a  shift  in  their 
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positions.  Thus  raanupulating  the  starting  signature  and  iterating, 
it  is  possible  to  get  back  the  specified  pressure  signature  at 
r/ l  =  0.5.  After  one  or  two  careful  iterations,  it  is  possible 
to  determine  the  body  of  revolution  for  the  types  of  specified 
pressure  signature  described  here. 

In  the  method  of  collocation  ±' (t)  is  assumed  to  be  of  the 


form 


f  (t)  = 


N 

V 

L . 
n=l 


ant 


The  choice  of  N  is  important ,  as  this  detem  ines  the  size  of  the 
matrix  ].  When  N  is  about  12,  there  is  loss  of  accuracy  in 
computation  because  of  the  type  of  matrices  occur ing  here.  To 
get  around  this  difficulty  and  to  use  the  "method  of  collocation" 
more  powerfully,  the  series  for  f ’  (t)  is  assumed  to  be  different 
in  different  portions  of  the  pressure  signature  and  collocation 
is  done  at  fewer  points  in  each  portion.  The  coefficients  of  the 
different  series  are  found  successively  in  the  calculation 
procedure.  Since  we  are  seeking  solutions  for  smooth  bodies, 
the  slopes  of  the  area  distribution  are  matched  at  the  junction 
points  in  the  calculation  of  the  coefficients  of  the  series. 

We  know  from  slender  body  theory  that 

f'(t)=s^i2 

2  IT 

Two  integrations  of  the  series  for  ff(t)  with  boundary 
conditions  S(o)  =  S*  (o)  =0  for  bodies  pointed  at  nose  give  the 
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area  distribution  and  thus  the  distribution  of  the  radius. 
Conclusions 

The  designed  bodies  of  revolution  which  produce  the  specified 
pressure  signature  at  Mach  number  M=1.2  and  a  spacing  r/l  =  0.5 
away  from  the  body  axis  have  been  plotted  in  Figure  2.2.9.  It  is 
seen  from  the  plots  that  the  body  (1)  of  revolution  due  to  the 
finite  rise  pressure  signature  is  much  sharper  at  the  nose  and 
thinner  than  body  (2)  due  to  the  plateau  pressure  signature.  This 
observation  could  otherwise  have  been  made  from  the  physical  point 
of  view,  but  the  present  analysis  gives  the  exact  shapes  cf  the 
bodies.  Although  in  the  present  example  problem  the  far-field  and 
near-field  theories  tend  to  merge  at  a  field  r/-t  =  4.2  and  M  =  1.2, 
for  thicker  bodies  these  cwo  theories  will  give  the  same  pressure 
signatures  at  distances  more  than  r/\{,  =  4.2  It  might  be  mentioned 
that  the  effect  of  Mach  number  on  both  theories  will  nearly  be  the 
same.  In  the  present  example,  it  has  been  checked  by  the  author 
that  although  the  pressure  signatures  due  to  the  two  bodies  de¬ 
generate  to  N-waves  at  large  distances,  body  (1)  produces  stronger 
front  shock  than  body  (2 ) . 

Recommendat ion 

Since  the  near-field  signature  shape  is  sensitive  to  slight 
deviations  in  the  shape  of  the  generating  body,  an-. extensive  study 

and  suitable  alterations  of  the  near-field  signature  will  make  it 
possible  to  reduce  sonic  boom  at  far-field  points  to  a  desired  level . 

It  also  may  be  pointed  out  here  that  the  same  theory  is  appli 
cable  for  a  body  of  revolution  at  an  angle  of  attack.  Although 
the  contribution  of  the  thickness  of  the  body  to  the  pressure 


signature  is  much  more  significant  compared  to  that 
small  angle  of  attack,  the  work  done  in  Secs.  2.2.1 
may  be  applied  to  determine  a  body  of  revolution  at 
of  attack  for  a  desired  pressure  signature. 
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y 

h  =  J  In 
o 


N  (r)  -  (y-t)1/2 
|N(r)  +  (y-t)1/2 


N (R)  +  (y-t) 1/2 
NCR)  -  (y-t) 1/2 


f'(t)  dt 


where 

N  (r)  =  N(yft,23r)  =  (y-t + 2pr)1/2, 


N(R)  =  N{y,t,2pR(y)}  =  {  y-t  +  2f3R(y)}  1/2 


The  integral  (A.l)  aiay  be  written  as 


L  =  L1**L2  +  L3“L4  ’ 

where 

y 

Lj  =  /  tn  {  (y-t  +  2pr)1/2  -  (y-t)1/2>  f'(t)  dt, 
o 

y 

L2  =  f  (y-t  +  2pr)172  +  (y-t)172}  f'(t)  dt, 
o 

y 

1*3  =  f  m  ({y-t +  2pR(y)}1/2+  (y-t)1/2]  f'(t)  dt, 
o 

y 

L4  =  /  tn[{y-t  +2pR(y)}1/2  +  (y-t)172]  f'(t)  dt 

o 


(A.l) 


(A. 2) 


(A. 3) 


(A. 4a) 


(A. 4b) 


(A .  4c  > 


(A. 4d) 


Integrating  by  parts  we  obtain  from  Eqs.  (a. 4) 


L,=f  (y)  in(2$r)V  2  -  f  (o)  ln{  (y+2pr) 1  y1/2}  / - r/^-L— 


2  o  (y-t)^2ty-t+2pr)V2 
(A.  5a) 


L2  =  f  (y)  -tn(23r)1/2-  f  (o)  *n{  (y+2pr)  1/2+y1/2}  +  -  /— • - -^1  d* - — 

2  2  o  (y-t)1/2(y-t+2Pr)1/2 


I3  =f  (y)^n[23R(y))1/‘  -  f  (oHn  [{  y+2£R(y  )}  A'  2+y1/2] 

y 

f  1  f _ f(t)  dt _ 

2  '  (y-t)1/2{y-t+23R(y)}1/2, 


(A. 5b) 


(A. 5c) 


L4  =  f  (y)tn{23R(y)}1/2-  f  (o)4n[{  y+2pfUy)}  1/2-  y1/2] 


_  1  J  3f(t)  dt _ 

2  (y-t )  1/2{  y-t+2pu  (y )}  1/2 


(A -5d) 


Combining  the  expressions  of  (A. 5)  an.i  cancelling  terms 
yields 


L  -  to  f {  <y+8Br>V^  ^  2H  (y+ZBBty))1  *  -yV2}]  f(o) 

L{  SOT  (y+26R(y))l/2+y^2}J 


y-t  +  23R(y) _ 

y-t+2£R (y )}  1/2 (y-t ) 1/2 


Assuming  the  body  to  be  pointed 


-y-t  +.2P.r _ 1  f(t) 

.-t+2Pr)1/2(y-t)1/2J 


(A.6) 


f  (o )  =  -S  =  R(o)  R'(o)  =  0, 


2.120 


I 

l 

I 


APPENDIX  2.2.2 


The  kernel  function 


\(z,t)  =  -  l  j  _ _ l1/2dj _ 

2  o  (1~^1/2{  (z-t)^+2pr}3/2 

Substituting  4=  sin2©  and  simplifying,  we  have 


(B.l) 


L,(z,t)  = - 1— 

(2pr)1/2 


0  (1+n  sinz©) 


de _ 

2l7l2"7l72 


(28r) 


_ d£ 

.2_.2„ 


Z“t  o  fl+n2sin2e)  (l-tn^sin2© ) 1/2  ’  (B‘2) 


where 


Z— t  9 

■  77“  =  n  >o  for  z  ^  t 
20r 


The  two  integrals  of  the  second  member  of  (B.2)  can  Le  reduced 
to  standard  forms  in  complete  elliptic  integrals  of  the  first  and 
second  hind  by  introducing  a  variable  *  in  place  of  e  by  means 


of  the  relation 


sn2v  =  (1-Bi2)  sin2© 

A  2  9  ) 

1+n  sin'2© 


(B.3) 


which  reduces  them  to  the  following  forms 


ir/2 

f—  de 

J  5 - 5 — ’ 

(1+n  sin  ©) 


K(<) 

•=  k'  f  d>;  =  k'kOc) 


2.122 


L  M 


K(jc) 


r _ de _ 

0  (l+n2sin2e)  (l+n2sin20)1//2 


=  k' $  dnxdX  =  k' E(k) 


where 


2  1/2 


K  =  (-*—) 
1+n 


,  0  <  JC  <  1, 

z-t+2pr 


/  ,,  2,1/2  ,  20r  7 

K  =  (1“/C  )  =  ( - *= - ) 

z-t-*-2pr 

is  the  complementary  modulus. 
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APPENDIX  2.2.3 


1.  Evaluation  of 


I,  (y)  =  / - ^  dt  _ _ 

Q  (y-t+2pr)1/2(y-t)2/2 

Changing  the  variable  of  integration  by  means  of  the  relation 
£=  (y-t)1/2,  I1(y)  reduces  to 
vl/2 

/  /  *2  m 

h<y)  =  2  f  SXZLJ. 

0  (r;  +2pr)1/2  (C 

Now  since 


(y-r)m  =  {Y  -  (£2+£{3r)}  m 

«  Y"1  -  (”>  +  ...  +  (-)r(")Y“-r(e2+2pr)r 

i  f  \ ®  —  .  m 


+  (-)  (£  +2pr)“ 


where 


Y  ^  y  +  2pr 


and,  furthermore,  since 


Tn  =  /  (52+2Br)m'2  dC  = 


2m- 1 


+  (22X1)  2flr 


,1/2 


m 


-I 


AL 


Ui2+2pr)1/2 


in 


(y^2  +  Y1/2 


(2pr) 


1/2 


) 


we  have  the  result  given  in  the  text . 


2.  Evaluation  of 


2m+l 


Vy)  =  f 1 - 172 - 172 

2  (y-t+2pr)1/2(y-t)1/2 


We  change  the  variable  of  integration  t  by  substituting 
2 

t  =  y  sin  e,  yielding 


Io(y)  = 


m+1 


ir/2 


(y+2pr) 


/ 


2  o+2 


sin  0  de 


1/2  J0 


(C.2) 


where 
2 


k-,  = 


y+2pr 


o 

,  0  <  k^  <  1  for  y  ^  0. 


Introducing  the  elliptic  function  sn  x  to  transform  the  variable 
of  integration  by  means  of  the  relation 


sn  X  =  sine 


the  integral  on  the  right-hand  side  of 


ir/2 


2ra+2 


-/ 


sin2ro+20  dO 
(1-k2  sin20)1/2 


K  (k,  ) 
/  sn 
o 


(C.2)  may  be  written  as 
2m'i'1X  dX  (C.3) 


In  reference  (Byrd  and  Friedman,  1954)  it  is  shown  that 
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2m+2 


2m(l+k1)  A2a  -  (2m-l)A2m^2 
(2m+l)  k? 


m 


| 

2,3,4 .  | 


and  the  first  two  integrals  Ag  and  are  easily  evaluated  to 
give 

K(kx) 

A 2  =  /  sn2X  dX  =  \  [K(R]L >  -  E^)], 
o  kl 

K(kx) 

A4  =  f  sn4x  dX  =  [  (2+k^)  K(kx>  -  2  (1+k^)  E(k1>J, 


where  K  and  E  are  complete  elliptic  integrals  of  the  first  and 
second  kind. 

The  above  formulas  will  enable  us  to  evaluate  I2(y)  in 
closed  form. 


j 

i 
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2.3  SONIC  BOOM  REDUCTION  BY  FOCUSED  LASER  BEAM  TECHNIQUES 

by  Mr.  Ronald  Kohl 
2.3.1  The  Problem 

If  a  supersonic  aircraft  could  be  made  sufficiently  long 
and  properly  shaped  and  designed,  more  favorable  ground  pressure 
signatures  could  be  obtained.  In  particular,  if  the  aircraft’s 
equivalent  body  of  revolution  has  certain  shapes,  the  over¬ 
pressure  can  be  reduced  or  the  shock  waves  in  the  signature 
can  be  eliminated  and  replaced  with  finite  rates  of  pressure 
rise.  The  latter  modification  in  a  signature  of  given  over¬ 
pressure  reduces  the  power  in  the  signature  which  occurs  in 
the  frequencies  audible  to  the  human  ear  (Hilton  and  Newman, 

1966;  Sepler  and  Harel,  1965;  Kryter,  1965). 

To  obtain  such  pressure  signatures,  however,  the  aircraft 
lengths  required  for  an  aircraft  of  350,000  lbs,  flying  at  Mach 
number  2.7,  at  an  altitude  of  65,000  feet,  range  from  500  ft 
up  depending  on  the  effective  bodyshape  desired  (McLean,  Carlson, 
and  Hunton,  1966).  Still  longer  lengths  are  required  for  the 
planned  SST  cruising  weight  of  about  600,000  lb.  Aircraft 
structures  of  these  lengths  are  considered  too  long  for 
economical  reasons.  A  means  must  be  found  to  create  airflows 
similar  to  those  that  would  occur  about  the  longer  aircraft, 
but  which  keeps  a  structural  aircraft  length  near  the  300  ft. 
length  considered  for  the  present  SST.  (See  Figure  2.3-1). 

To  deflect  the  airflow  without  using  structure,  energy  could  be 
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aircraft 


introduced  into  the  airflow  and  here  this  energy  must  be  intro¬ 
duced  at  distances  the  order  of  100  ft  or  more  from ^the  aircraft 
structure.  To  do  this,  novel  means  must  be  used. 

The  idea  was  advanced  to  use  present  or  future  lasers  to 
accomplish  this.  The  feasibility  of  the  use  of  lasers  depends 
on,  among  other  things,  the  power  required  to  be  introduced  into 
the  airflow.  Of  particular  interest  is  the  power  required  in 
the  region  out  in  front  of  the  aircraft  where  the  front  of  the 
pressure  signature  would  be  affected  and  where  the  introduction 
of  burning  fuels  may  not  be  feasible. 


In  the  presence  of  one  of  these  uneconomica lly  long  noses 
of  a  desired  shape,  the  flow  streamtubes  take  certain  forms. 

The  flow  outside  of  a  streamtube,  including  the  ground  pressure 
signature,  is  completely  determined  by  the  streamtube.  No  matter 
|  what  means  are  used  to  produce  the  streamtube,  if  the  streamtube 

[  is  produced,  the  desired  ground  pressure  signature  is  obtained. 

The  idea  here  is  to  produce  the  streamtube  desired  by  adding  heat 
to  the  flow  inside  the  streamtube  rather  than  having  a  nose’  of 
the  desired  shape  located  inside  the  streamtube.  In  this  why  there 
is  a  "phantom"  nose  present.  To  find  the  power  requirements,  we 
find  the  streamtube  area  developments  and  the  pressures  produced 
on  the  surfaces  of  the  streamtubes  by  the  uneconomica lly  long 
nose  of  a  desired  shape.  The  input  power  required  to  duplicate 
the  area  development  and  surface  pressures  of  a  given  streamtube 
is  then  found  by  considering  the  streamtube  as  a  one  dimensional 
1  channel  and  finding  the  heat  addition  required  to  duplicate  the 

S'  surface  pressures  which  exist  in  the  actual  flow.  For  a  given 

nose  the  calculation  is  done  for  several  streamtubes  of  differing 
initial  area  to  determine  the  effect  on  the  streamtube  power  re¬ 
quirement  of  the  initial  streamtube  size. 

To  obtain  linear  finite  rates  of  ground  signature  pressure 
,  rise  with  a  structural  nose,  application  of  Whitham's  approach 

(Whitham,  1952)  shows  the  well  known  axisymmetric  area  development 
F  S  =  ax5^^  is  required  where  x  is  the  axial  distance  from  the  nose 

tip.  The  constant  a  is  related  to  the  desired  rate  of  ground 
signature  pressure  rise  dP/dt  by 


2. 130 


a 


16  vA3„r 


15M2  [!KIP  U  /(dP/dt)  +  (i'+l)  M  rp'1] 

where  M^,  and  are  the  upstream  Mach  number,  static  pressure, 
and  velocity  at  the  altitude  of  interest  r.  The  quantity  (3  is 

_  CO 

/2  ^ 

/Mm  “  1  and  y  is  the  ratio  of  the  specific  heats  for  air.  is 
the  reflection  coefficient,  the  factor  by  which  the  magnitude  of 
the  pressure  in  the  ground  pressure  signature  is  increased  over 
the  free  stream  pressure  signature  due  to  the  presence  of  the 
ground.  Here  K  -  2  was  used.  K  is  the  correction  coefficient 

^  3 

of  Kane  (Kane  1966;  Kane,  1967)  which  is  the  ratio  of  the  pressure 
variation  at  sea  level  in  the  U.  S.  Standard  Atmosphere  to  the 
pressure  variation  at  sea  level  in  a  uniform  atmosphere  of  pressure 

P  .  For  an  altitude  of  60,000  ft.  this  factor  is  K  «  4.33.  For 

oo  a 

the  purpose  of  this  work  the  small  differences  in  pressure  signature 
lengths  between  the  uniform  and  standard  atmosphere  signatures  are 
ignored.  (Hayes  and  uaefeli,  1968). 

To  find  the  streamtube  area  or  radius  development  with  x 
and  the  streamtube  surface  pressures  for  a  given  nose,  the  nose 
was  replaced  by  flow  sources  on  its  axis  in  the  usual  slender 
body  -  perturbed  flow  approach.  The  stream!' unction  Y,  when  written 
in  t.he  form 

*'pa,  =  V +  *> 

where  is  the  upstream  density  at  the  altitude  of  interest  and 
r  is  the  radial  distance  from  the  axis,  was  found  to  be  given  by 
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*  = 


-1 

cosh 
-  / 


f (x-pr  cosh  c ) 


pr  coshn  d< 


o 

where  f(x)  is  the  usual  axial  source  distribution  in  the  slender 
body  approach 


f  (x) 


_  1_  dS 
2ir  dx° 


The  radius  of  a  streamtube,  selected  upstream  where  f  =  0,  was 
then  traced  through  the  flow  using  the  constancy  of  Pressures 
on  the  streamtube  surface  were  found  through  the  pressure  coefficient 


2u  -  v2 , 


where  the  flow  velocity  in  the  axial  direction  is  (1  +  u)  Uo,  and 

in  the  radial ,  vU  ,  with  u  and  v  given  by  the  well  known  relations 

00 

cosh_1(|-) 

u  =  -  J  f’(x-pr  cosh o)  do 

o 

cosh-1  (^p) 

v  =  J  f’(x~Pr  cosbo)  p  coshc  do. 

o 

where  f ’ (x)  =  df/dx.  In  order  to  calculate  the  heat  addition  re¬ 
quired  to  obtain  pressures  in  the  streamtube-channel  which  match 
the  known  pressures  at  the  surface  of  the  streamtube,  the  inviscid 
equations  of  conservation  of  mass,  momentum  and  energy,  the  equation 
of  state  (perfect  gas)  and  the  definition  of  the  stagnation  tempera¬ 
ture  T  are  combined  under  the  assumption  of  one  dimensional  depen- 
o 

dence  to  obtain  the  power  addition  per  unit  axial  length 
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dP 

dx 


p  U  A 

rco  00  00 


where 


=  [1  +  <*-!>  M2/2]_1  jr  jj§ 


+  p2  ir -1m-2 [  1  +  (r-i)  m2/2]"1  i  g| 

2  2  2  2 
with  p  =  M  -1  and  M  the  square  of  the  Mach  number.  M  is  found 

down  the  channel  by  applying 

15T  =  -P‘2f1+  <»-»  m2/2]  <1  +>“2>  F* sir 

M  O 

+  2(3_2[l  +  or- 1)  112/2]  i  H 

These  expressions  can  be  found  in  Table  8.2  of  (Shapiro,  1953). 

The  calculation  moved  down  the  channel  by  steps  of  axial 
displacement  ax,  the  approximation  to  dx,  which  was  chosen  to  be 
sufficiently  small  for  the  sake  of  accuracy,  but  not  so  small  as 
to  take  unreasonable  computer  time.  The  quantity  dA/A  was  used  in 
the  form  2dr/r  where  dr  is  known  from 


i 


i 


■  1  CO 


where 
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M 

or 


-/ 


COSh 

Pr  /  2  2 

f  (x-pr  cosh  o )  p  r  cosh  a  d.a 


-  J 


cosh-1  (~~) 


f  (x-pr  cosh  a)  p  cosh  o  da 


and 


dx 


-/ 


cosh'1^) 

f '  (x-pr  cosh  a  )  pr  cosh  a  da . 


The  change  in  streamtube  radius  Ar  with  axial  displacement  ax 
from  axial  position  x^  was  found  by  using  Ar  from  a  previous 
trial  (beginning  with  the  final  value  of  Ar  from  the  previous 
step)  to  obtain  ra  =  ri  +  Ar/2.  This  value  was  used  with 
xa  =  x-^  +  ax/2  in  cty/dr  and  <fy/dx  to  obtain  a  new  trial  Ar. 

This  process  was  continued  until  the  relative  difference  in 

g 

succeeding  trial  Arfs  was  less  than  1:10  .  The  streamtube 
radius  at  x2  =  xi  +  Ax  was  then  rg  =  r^  +  Ar.  To  guard 
against  wandering  from  the  original  streamtube  in  moving  down 
the  flow,  the  difference  in  the  value  of  the  actual  and  original 
stream  function  was  monitered.  The  relative  error  in  stream- 

3 

tube  radius  was  always  less  than  1:10  . 

The  quantity  dp/p  is 


i 


Where  the  nose  first  begins  to  influence  the  streamtube  aCd 
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is  the  order  of  Cp  so  Ap/p,  wfcr.ch  was  alv.ays  small ,  was  used  in  the 
form 


A£  =  „ 

P  c 


2  (C 


S2 


Aieii 


P2  +  cPi  + 


4/7  M‘ 


In  moving  down  the  stre2mtube,  the  numerator  was  monitered 
for  possible  loss  of  significant  figures  as  Cpj*  approached 


Knowing  aA/A  and  Ap/p,  ATq/T0  can  be  determined  accurately 
by  knowing  M2  at  x=xa>  For  small  ax,  Ma  =  0^  +  Mg)/2  and  Mg 
is  found  by  successive  trials,  beginning  with  a  trial  value 
for  M  (initially  M-.  was  used),  obtaining  a  trial  AT  /T  and 

a  1  O  O 

2  2 

AM  /M“  and  thus  a  new  value  for  Mg  and  M&.  This  was  continued 

S 

until  successive  trial  values  of  M„  differed  by  less  than  1:10  . 

a 

This  approach  in  obtaining  M  was  also  used  in  (Miller  and 

A 

Carlson,  1969). 

A  program  was  developed  to  calculate  the  desired  quantities 
along  the  lines  given  above  for  phantom  noses  of  area  develop¬ 
ment  S  =  ax11  with  n^2.  All  integrands  for  such  noses  are 
finire.  The  integral  evaluations  in  this  program  were  done 
with  a  combination  of  Simpson’s  rule  and  Newton’s  3/8  rule 
with  the  number  of  intervals  adjusted  for  accuracy  and  speed. 


To  speed  up  the  computer  runs  on  conical  noses  h  =  2)  the 

integrals  were  done  analytically.  The  results  obtained 

using  numeric  and  analytic  integrations  always  agreed  to  better 
3 

than  1*10  .  All  computer  runs  were  don^  on  the  IBM  1130 
system  here  at  the  Institute. 

To  establish  minimum  power  requirements,  noses  which  pro¬ 
duce  finite  pressure  rises  in  the  pressure  signature  were 
joined  tangentially  to  the  effective  body  of  revolution  of  an 
SST  so  that  the  body  was  enclosed  within  the  nose.  See  Figure  2.3 -2. 
The  power  required  to  duplicate  the  flow  caused  by  the  presence 
of  the  phantom  nose  between  the  phantom  nose  tip  and  the  SST 
body  tip  was  then  calculated  by  summing  the  power  require¬ 
ments  per  step  down  the  channel  until  the  characteristic  from 
the  SST  body  tip  was  reached.  To  duplicate  the  flow  produced 
by  the  phantom  nose  behind  the  position  of  the  SST  body  tip  charac¬ 
teristic  would  require  cooling  the  flow  to  compensate  for  the  com¬ 
pression  of  the  flow  impinging  on  the  SST.  If  no  cooling  is 
allowed  and  we  still  wish  flow  behind  the  SST  b  tip  to  have 
the  characteristics  of  flow  over  the  phantom  nose,  an  even 
longer  phantom  nose  would  be  required  for  a  given  rate  of  pres¬ 
sure  signature  rise.  Thus  for  a  given  nose  shape,  the  length 
of  nose  considered  here,  the  distance  from  the  phantom  nose 
tip  to  the  SST  body  tip,  is  a  minimum  and  the  power  requirements 
obtained  are  minimum  power  requirements. 
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2.3.3  Results 


Results  for  a  typie.- ;  case  are  given  in  Figures  2.3.  3, .4,.  5 

and  .6.  In  this  case  a  phantom  nose  which  produces  a  ground 

2 

pressure  rise  of  2  lb/ft  in  10  msec  from  an  altitude  of  60,000 

feet  at  a  Mach  number  of  2.7  is  joined  to  the  effective  body 

of  a  600,000  lb.  SST  resulting  in  a  phantom  nose  length  of  249 

feet  to  be  duplicated  by  heat  addition.  In  this  example  the 

2 

streamtube  of  n00  ft  initial  area  is  the  one  considered  as 
shown  in  Figure  2.3-3.  In  Figure  2.3.4  the  required  power  in¬ 
put  per  unit  axial  length,  the  required  input  power  distribution, 
is  plotted  down  the  streamtube.  Figure  2.3-5  shov/s  the  integral 
of  the  curve  in  Figure  2.3-4,  or  the  total  power  required. 

Using  the  energy  available  from  the  combustion  of  a  pound  of 
jet  kerosene  and  the  SST  cruise  velocity  of  1,780  m/hr,  this 
quantity  is  given  in  pounds  of  fuel  per  mile  of  sonic  boom 
abatement  system  operation.  The  power  required  to  simulate 
the  presence  of  the  249  foot  nose  is  the  power  produced  by 
combusting  43  pound.s  of  fuel  per  mile  or  420  Megawatts.  Figure 
2.3.6  shows  that  this  result  depends  only  slightly  on  the  stream- 
tube-channel  chosen  for  the  cross  sectional  heat  addition.  The 
input  power  distribution  is  also  practically  independent  of 
streamtube  area  as  shown  in  Figure  2.3.7. 

A  comparison  of  the  power  requirements  for  various  phantom 
noses  is  shown  in  Figure  2.3.8.  The  power  requirements  are 
plotted  against  the  rate  of  ground  signature  pressure  rise  pro¬ 
duced  by  the  noses.  These  rates  range  from  a  barely  finite  rate 
to  p  rate  where  the  ground  pressure  does  not  reach  2  lbs/ft2,  the 
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unmodified,  predicted  SST  maximum  overpressure,  for  100  msec, 
which  is  one  fourth  the  predicted  SST  pressure  signature  length. 
Thus,  Figure  2.3.8  indicates  that  if  a  barely  finite  rate  of 
pressure  rise  can  be  obtained,  beneficial  rates  of  rise  can  be 
obtained  with  an  increase  in  power  of  the  order  of  50%  or  less. 

A  comparison  of  power  requirements  as  a  function  of  nose 
shape  was  also  made  to  see  if  relaxing  the  phantom  nose  shape 
from  the  desirable  S  =*  ax5/"2  shape  might  bring  a  helpful  re- 

2 

duction  in  the  power  requirement.  Here  a  conical  nose  (S-ax  ) 

3 

and  a  nose  of  area  development  S  =  ax  ,  both  of  the  same  length 
and  base  area  as  the  nose  of  Figure  2.3.3  are  compared  to  the 
nose  of  that  figure.  A  longitudinal  cross  section  of  the  noses 
is  shown  in  Figure  2.3.9.  The  resultant  input  power  distributions 
are  compared  in  Figure  2.3.10  for  an  initial  streamtube  area  of 
100  square  feet.  A  comparison  of  the  total  power  required  tc 
duplicate  the  effects  of  these  noses,  and  noses  similarly  ob¬ 
tained,  is  shown  in  Figure  2.3.11.  While  the  input  power  dis¬ 
tributions  are  some  what  different,  the  total  power  required 
is  practically  the  sam.  for  the  various  shapes.  (One  can  show 

that  for  slender  noses  of  area  development  S  =  axn,  the  linear 
5/2 

rise  or  S  =  ax  nose  shape  produces  a  minimum  average  pressure 

rise  and  a  minimum  maximum  pressure  rise  in  a  shockless  signature, 

as  well  as  being  the  shape  which  produces  a  shockless  pressure 

signature  at  a  maximum  distance  from  the  nose  axis-  A  shape  of 

S  =  ax11  with  n  less  than  5/2  does  not  produce  a  shockJess  signature). 

2  3 

The  power  requirements  for  the  S = ax  and  S = ax  noses,  as  with 
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10  FT  LB.  SEC 


FT. 


FI(  URE  2.3.4.  Input  power  per  unit  axial  length.  The 
input  power  distribution  required  to  duplicate  the 
nose  and  strearatube  oi  Figure  2.3.3. 
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SST  CRUISE  REQUIREMENT 


LBS.  QF  p] 


"ITIaL  STSEAlfTUBE  ARM 
F«iURE  2.3.6.  powe„  „  . 

°r  f » reamt ube  i“rtSt|,VS;  initial 

"  FRantom  nose  of  Mg„|V?a|  ls  added 


LBS.  OF  FUEL 


MILE 


RATE  OF  PRESSURE  RISE  PRODUCED  LE/FT2 

BY  PHANTOM  NOSE  msec 


FIGURE  2.3.8.  Power  requirement  vs.  rate  of  ground 
signature  pressure  rise  produced  by  phantom 
nose  for  SST  of  600,000  lbs  at  Mach  number  2.7. 
(The  length  of  nose  duplicated  is  indicated  in 
parentheses  in  feet.) 
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(LINEAR  RISE) 


S  ax  NOSE 


ax2NOSE 

(CONE) 


250  FT. 


FIGURE  2.3.9.  Comparison  of  nose  shapes.  Longitudinal 
cross  sections  of  noses  ol  the  same  length  and  base 
area.  The  linear  rise  nose  is  the  nose  of  Figure 
2.3.3. 
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1 06  FT . 


LD./Siv 


F' 


/ 

/ 


/ 

/ 


AXIAL  DISTANCE  FROM  PHANTOM  NOSE  TIP 


FIGURE  2.3.10  A  comparison  of  the  input  power  distri¬ 
butions  required  to  duplicate  the  linear  rise  nose 
of  Figures  2.3.3  through  2.3.6  and  the  (S  -  ax3/2) 


and 


S  -  ax^  noses  of  the  same  length  and  base  area. 
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DBS.  OF  FUEL 
MILE 


LINEAR  RISE  NOSES 


S  -  ax  NOSES 
OF  SAME  LENGTH 
AND  BASE  AREA 


CONICAL  NOSES  OF  SAME 
LENGTH  AND  BASE  AREA 


INITIAL 
STREAMTUBE  r 
AREA  =  100  FT* 


RATE  OF  PRESSURE  RISE  PRODUCED 
BY  LINEAR  LINE  NOSE 


LB/ FT 
mSEC 


FIGURE  2.3.11.  A  comparison  of  the  power  required  to 
duplicate  the  linear  rise  (S  -  ax5/2)  noses  of 
Figure  2.3.8,  and  the  conical  (S  -=  ax2)  and  S  -  ax3 
noses  of  the  same  length  and  base  area. 
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MACH 

NUMBER 


FIGURE  2.3.13=  Streamtube  -  channel  Mach  number  for  the 
phantom  ncses  of  Figure  2.3.9. 
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the  linear  rise  noses,  do  net  depend  very  much  on  the  particular 
streamtube-channel  chosen  as  shown  in  Figure  2.3.12.  For  all 
the  nose-st reamt ube  combinations  the  Mach  number  in  the  stream- 
tube-channel  stayed  above  one.  See  Figure  2.3.13. 

Consideration  of  the  above  results  leads  to  an  approxi¬ 
mation  scheme.  In  the  expression  for  TQ^  dTQ/dx  since  the 
p  *dp/dx  term  is  observed  to  be  small  compared  to  the  A  ^dA/dx 
term  except  near  the  region  where  the  streamtube  considered  is 
first  affected  by  the  nose,  and  as  M  and  TQ  undergo  small  changes 
in  that  region  compared  to  the  changes  occuring  further  down  the 
streamtube,  we  neglect  the  effect  of  the  p  1  dp/dx  term  and  write 

dT 

— ~  Sr  [1  +  (•y-'l)M2/2]"1  — 

T  A 

o 

and 

dM2  ~  _  dA 
M2  A 

2 

Thus  M  A  is  a  constant  which  can  be  evaluated  upstream  to  give 
M2  =  M“  A  /A  where  A  is  the  upstream,  or  initial,  streamtube 

oo  oo  co 

area .  Then 


dT 


-  =  [A  +  (7-1)  M2A  /2]"1  dA 


,-l 


and  T  fA  +  (-v-l)  W“A  /2l  is  a  constant.  This  uplfes  that  the 

Ol  '  w  00  J 

power  addition  per  unit  streamtube-cha»inel  length  is 
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dP  _  ^oo^oo*'  Oa>  dA  _  u  c  T  ^ 
dx  l+(7-l)M2/2  dx  P°°  00  p  00  dx 

where  is  the  upstream  static  temperature.  Since  the  stream- 
tubes  considered  here  are  sufficiently  close  to  the  axis  com¬ 
pared  to  the  phantom  nose  surface,  slender  body  theory  should 
apply  and  the  rate  of  streamtube  area  change  at  a  given  cross 
sectional  plane  should  be  determined  to  a  sufficient  approxi- 
aation  by  the  axial  flow  sources  at  the  intersection  of  tne  ■ 
plane  and  axis.  That  is, 

—  =  2ir  f  (x)  =  ~ 

dx  dx 

where,  as  above,  S  is  the  cross  sectional 
nose.  The  channel  heat  addition  per  unit 
written  as 

—  =  ,U  c  T  — . 

.  1  oo  oo  p  CO  . 

dx  ^  dx 

This  last  expression  essentially  reproduces  the  curves  of 
figure  2.3.10  except  for  the  spike  at  the  front  of  the  conical 
nose  curve  which  is  due  to  the  steep,  but  small,  increase  in 
pressure  there,  a  feature  one  would  not  wish  to  duplicate  in 
any  event.  This  last  expression  was  derived  by  Sv.igart  and 
Lubard  (Swigart  and  Lubard,  1969)  in  a  different  approach  and 
the  above  discussion  shows  the  circumstances  under  which  the 
answers  obtained  by  their  approach  will  agree  with  those  obtained 


area  of  the  phantom 
length  can  then  be 
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by  a  one  dimensional  channel  approach  as  used  here  and  else¬ 
where  (Miller  and  Carlson,  .1939;  Siegelman,  1970).  Relating 
the  mass  source  per  unit  length  represented  by  the  axial  sources, 
S,  =  2ir o  U  f  (x)  ,  to  the  equivalent  heat  sources  per  unit 

1  *  m 

length.  Q.  =  c  T  S.  (Tsien  and  Beilock,  1949;  Wiilmarth,  1957), 

i  p  ®  i, 

one  has  dP/dx  =  Q  to  the  above  approximation. 


2.3.4  Discussion 


The  power  requirements  as  obtained  above  have  serious  con¬ 
sequences  for  the  use  of  any  laser  system.  For  the  purpose  of 
this  discussion  we  shall  use  rhe  420  Megawatt  power  require¬ 
ment  figure  for  the  phantom  nose  of  Figure  2.3.3  which  gives  a 

2 

finite  rate  of  pressure  rise  of  0.2  (lb/ft  )/m  sec,  equivalent 

o 

to  the  attainment  of  2  lb/ft*"  in  10  milliseconds.  The  power 
requirement  for  this  rate  of  rise  is  compared  to  those  of  other 
rates  of  rise  in  Figure  2.3.8. 

Many  laser  systems  rely  on  electric  energy  as  their  initial 
energy  source  either  through  the  use  of  flashlamps  or  through 
electrical  discharges  of  various  types.  Questions  of  laser 
efliciency  aside,  the  generation  cf  420  Megawatts  of  electrical 
power  carries  a  severe  weight  penalty.  At  these  power  levels 
magnetohydrodynamic  (MriD)  generators  offer  the  smallest  equipment 
weight  per  kilowatt,  but  turbo-alternators  offer  the  lowest 
specific  fuel  consumption  assuming  environmental  air  is  used  as 
the  oxydizer  (Cooper,  1971).  The  '!HD  equipment  weight  is  about 
0.2  pounds  per  kilowatt  or  84  thousand  pounds  for  the  generation 
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of  420  Megawatts.  This  weight  is  30%  of  the  nonfuel  takeoff 
weight  of  the  commerical  version  of  the  Boeing  SST  (Aviation 
Week,  1970).  The  MHD  geneiator  specific  fuel  consumption  is 
about  1. 1  (lb/sec)/MW  or  930  lb  of  fuel  per  mile  of  sonic  boom 
abatement  operation  during  SST  cruise.  With  the  turbo-alternator, 
which  Will  involve  higher  equipment  weight,  this  figure  is  470 
lb  of  fuel  per  mile.  Both  figures  assume  i he  use  of  environmental 
air  as  the  oxydizer. 

There  are  laser  systems  recently  announced  in  the  literature, 
which  do  not  require  the  generation  of  electric  power.  These  are 
the  purely  chemical  or  direct  combustion  laser  (Cool  and  Stephens, 
1969)  and  gas  dynamic  laser  (Gerry,  1970;  Meinzer,  1971).  Putting 
questions  of  suitability  for  heating  the  airflow  aside,  these 
carry  severe  weight  penalties  also.  The  chemical  laser  referenced 
here  requires  no  electric  discharge  to  attain  partial  dissociation 
of  the  reactants,  a  procedure  which  is  used  on  other  chemical 
lasers  and  typically  requires  an  electrical  input  power  many  times 
the  laser  output  power  (Spencer  et.  al. ,  1970).  But  the  laser 
type  referenced  here  do*'?  require  a  non-cycling  mass  flow.  Accord¬ 
ing  to  the  article  cited,  this  How  would  have  to  be  so  large  that 
any  practical  applxcatian  is  fully  unfeasible.  There  are  unclassi¬ 
fied,  unconiirraed  reports  of  the  attainment  or  future  attainment 
of  100  kilowatts  per  pound  of  flow  per  second  from  a  direct  com¬ 
bustion  laser,  but  this  would  still  leave  an  unfeasibly  large 
flow  requirement. 

With  gas  dynamic  lasers  the  spscific  power  figure  of  2  to  3 
kilowktts  per  pound  of  flow  per  second  is  still  far  too  small 

for  a  single  use  flow.  In  gas  dynamic  lasers  the  mass  flow  can 
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be  cycled.  The  laser  output  would  then  originate  in  the  com¬ 
pressor  driving  the  flow  circuit  (Hertzberg,  Johnston,  and 
Ablstrom,  L9711.  The  best  mechanical  power  generation  to  equip¬ 
ment  weight  ratio  has  recently  movec  from  5.0  kW/lb  (Wood,  1968) 
to  6.5  kW/lb  (Aerospace  Association,  1970).  Even  at  the  latter 
figure,  however,  65,000  pounds  of  equipment  24%  of  the  nonfuel  SST 
takeoff  weight,  would  be  required  to  generate  420  Megawatts  of 
mechanical  power  for  the  compressor.  This  figure  does  not  include 
the  weight  of  the  compressor  itself.  It  is  to  be  kept  in  mind 
that  the  sonic  boom  reduction  by  means  of  lasers  can  be  readily 
restricted  to  sensitive  parts  of  the  overflight  terrain  by  means 
of  turning  on  and  off  the  laser  equipment. 


2.3.5.  Conclusions  and  Recommendations 

It  is  concluded  that  large  amounts  of  power  will  be  required 
to  produce  desirable  finite  rise  time  pressure  signatures.  This 
is  shown  m  subsection  2.3.3  It  is  also  concluded,  as  discussed 
in  subsection  2.3.4,  that  at  present  there  exists  no  system  which 
can  generate  such  power  in  laser  emission  without  severe  weight 
penalties. 

In  view  of  the  amazing  advance  in  the  maximum  average  output 
power  of  lasers  in  this  past  decade,  laser  progress  should  be  infor¬ 
mally  monitored  for  the  discovery  of  a  system  which  obtains  its 
output  chemically  from  its  fuel  with  power  output  to  mass  flow 
ratios  approaching  those  available  from  the  combustion  of  commercial 
fuels.  It  should  be  borne  in  mind  that  heating  of  the  airstream 
by  laser  emission  can  be  accomplished  not  only  by  gas  breakdown, 


2.  4.  ENGINE- AIR FLAME  INTEGRATION  WITH  SPECIAL  EMPHASIS 

ON  NON-CIRCULAR  ENGINE  EXHAUSTS  AND  JET  FLAPS 

2.  4.  1  Pressure  Distribution  of  Deflected  Two-Dimensional  Jets  Behind  Wings 

by  Dr.  B.  H.  Goether';,  assisted  by  Dr.  Heinz  Gruschka  and  UTSI 
Student  Mr.  Philip  Kessel 

One  major  part  of  the  sonic  boom  signature  is  produced  by  the  lift  of 
an  aircraft;  this  part  is  usually  much  larger  than  the  sonic  boom  part  due  to 
the  solid  displacement  of  the  aircraft  structure.  Thus,  the  major  effort  of  this 
study  was  directed  towards  determining  which  sonic  boom  alleviations  can  be 
obtained  by  appropriately  utilizing  the  lift  component  of  deflected  exhaust  jets 
of  turbojet  or  ramjet  engines. 

On  Figures  2.4.  1  to  2.4.4  some  fundamental  relationships  in  the  forma¬ 
tion  of  the  sonic  boom  signature  for  two-dimensional  lifting  surfaces  are 
depicted.  It  is  demonstrated  on  Figure  2.4.  1  that  the  downward  momentum 
behind  a  wing  without  ground  effect  corresponds  to  one-half  of  the  total  lift  on 
either  side  of  the  wing.  On  the  other  hand,  if  the  pressure  and  expansion  waves 
impinge  upon  the  ground,  the  total  impulse  submitted  to  the  ground  is  equal  to 
the  full  amount  of  the  lift.  As  Figure  2.4.2  also  indicates,  the  dividing  stream¬ 
line  behind  the  wing  approaches  again  the  undisturbed  position  after  passing  through 
the  reflected  wave  system,  because  all  waves  cancel  each  other  at  large  distances. 

On  Figures  2.  4.  3  and  2.  4.  4  it  is  demonstrated  that  the  sonic  boom  signa¬ 
ture  on  the  ground  spreads  over  a  larger  distance  and  has  a  smaller  pressure 
peak  when  the  lifting  wing  has  a  larger  chord;  that  is  when  the  lift  is  spread 
over  a  larger  distance  in  the  flight  direction. 
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The  purpose  of  this  investigation  is  to  show  that  the  distribution  of  the 
lift  in  the  flight  direction  can  occur  not  only  by  means  of  stretching  the  wing 
chord,  but  also  by  having  a  part  of  the  total  lift  produced  by  the  deflected  exhaust 
jet  of  the  engine  (see  Figure  2.4.  5).  In  the  case  of  a  deflected  exhaust  jet, 
it  is  to  be  understood  that  the  entire  vertical  component  of  the  deflected  jet, 
that  is  the  lift  component,  is  transferred  to  the  wing  at  the  deflection  point 
of  the  jet,  that  is  in  the  exhaust  nos?  ;<•.  However,  as  far  as  the  action  on  the 
flow  and  particularly  also  on  the  wave  system  which  finally  impinges  on  the 
ground  it.  concerned,  the  lift  component  of  the  jet  is  distributed  over  a  distance 
downstream  of  the  exhaust  nozzle. 

This  phenomenon  has  not  been  fully  recognized  in  the  past.  Therefore, 
on  Figure  2.4.6,  various  force  equilibrium  sketches  are  prepared  which 
show  that  the  reaction  on  the  flow  behind  the  wing,  equivalent  to  the  vertical 
lift  component  of  the  jet,  is  caused  by  the  aerodynamic  pressure  difference 
between  the  upper  and  the  lower  surfaces  of  the  deflected  jet.  In  effect,  the 
deflected  jet  acts  like  a  solid  plate  in  the  flow  which  is  curved  according  to  the 
jet. 

A  theoretical  two-dirr  ional  calculation  was  conducted  as  shown  in 
principle  on  Figure  2.4.7  and  on  pages  2.  167  a  and  b. 

A  numerical  evaluation  of  the  pressure  distribution  curves  is  shown  on 
Figure  2.4.  8.  The  rearward  shifting  of  the  lift  is  larger  the  larger  the  Mach 
number  is  and  the  larger  the  jet  momentum  coefficient  is. 
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A  comparison  of  the  jet  lift  was  conducted  with  respect  to  the  lift 

generation  by  means  of  wings  and  mechanical  flaps  (see  Figure  2.4.9). 

A  jet  flap  operating  at  a  Mach  number  of  M  =  3  with  a  momentum  coefficient 

of  C  =1.4  has  the  same  lift- producing  effectiveness  as  a  wing  in  supersonic 
M- 

flow.  In  subsonic  flow,  a  jet  momentum  coefficient  of  approximately 

C  =  1.  9  is  required  to  produce  the  same  lift  as  a  wing. 

M- 

At  a  jet  momentum  coefficient  of  C  =  .  5,  figure  2.4.9  shows  that 

M- 

in  flight  at  Mach  number  3,  the  jet  flap  would  be  identical  In  effectiveness 
as  a  mechanical  flap  with  a  flap  chord  of  approximately  35  percent  of  the 
main  wing.  The  same  jet  momentum  coefficient  C  =  .  5  in  subsonic  flow 
would  have  the  same  lift  effectiveness  as  a  flap  of  13  percent  chord  length. 

In  the  following,  some  applications  of  the  new  theory  on  lift  shift 
and  performance  changes  due  to  jet  flap  deflection  are  shown  for  the  Mach 
number  range  up  to  M  =  12. 


JET  FLAP  IN  SUPERSONIC  FLOW 


JET  FLAP  IN  SUPERSONIC  FLOW 


2.4.2  SHIFT  OF  LIFT  BEHIND  SOLID  WING  AND  AIR CRA  FT  PERFORMANCE 
CHANGES  DUE  TO  DEFLECTED  JETS 


Figure  2.4.  10  shows  the  airframe  and  propulsion  system  data  used  in 

this  study.  The  Lift  Over  Drag  (L/D)  ratios  from  various  configuration  studies 

were  plotted  and  a  boundary  curve  drawn  which  is  believed  to  be  representative 

of  the  most  optimistic  L/D  ratio  obtainable.  The  curves  in  Figure  2.4.  10  give  also 

ramjet  and  turbojet  performance  and  the  values  of  air  specific  impulse,  I  .  , 

as  a  function  of  Mach  number.  This  data  is  used  to  calculate  the  exhaust 

momentum  to  thrust  ratio  and  also  the  jet  momentum  coefficient,  C  . 
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J  =  exhaust  momentum 

m 

v  =  flight  velocity 

CO 

g  =  gravitation  constant 

T  =  thrust 

I  .  =  air  specific  impulse, 

air 

The  data  from  Figure  2.4.  10  was  used  to  calculate  the  ratio  of  jet  lift  to  wing 
lift.  This  data  is  plotted  in  Figure  2.4.  3  1.  The  effect  of  jet  flap  on  aircraft 
L/D  and  range  is  shown  in  Figure  2.4.  12.  In  these  plots  the  advantage  of  exhaust 

deflection  at  high  Mach  numbers  shows  quite  plainly.  It  is  also  obvious  that  for 
every  Mach  number  there  is  an  optimum  exhaust  deflection,  measured  from  the 
horizontal,  which  will  maximize  range.  The  optimum  range  factor  varies  between 
+0.  7%  at  6°  for  Mach  2,  to  +27%  at  13.  5°  for  Mach  12. 
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The  expression  for  pressure  distribution  along  a  jet  .flap  can  be  rearranged 
in  terms  of  aircraft  and  propulsion  parameters  by  substituting  the  engine  mass 
flow  and  aircraft  geometry  factors  we  get: 


B.*P  _  B.  r  4q_  J 

_LX  -  J-  .  W/8 

te  S  te 


e  L  a 
te 


mu  te 
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where 


BaP 


B  W/S 
te 


Non-dimensional  flap  lift 


Width  of  jet  flap 


Trailing  edge  span 
aircraft  wing  loading 
pressure  across  flap  at  status  x 
jet  flap  deflection 
dynamic  pressure 

Vm 

CO 

engine  mass  flow 
engine  exhaust  velocity 


The  pressure  distribution  for  some  typical  aircraft  are  shown  in  Figure  2.4.  13, 


2.  4.  14,  and  2.  4.  1 5. 

A  lift  discontinuity  exists  at  the  trailing  edge  of  the  wing  due  to  the  abrupt 
change  in  the  direction  of  the  airflow  caused  by  ine  jet  flap.  The  extent  of  this 
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discontinuity  can  I.e  limited  by  decreasing  the  span  of  the  jet  flap.  Since  momentum 
considerations  insure  that  the  total  lift  must  remain  the  same,  we  can  predict  that 
the  lift  distribution  of  the  jet  with  smaller  span  must  extend  further  behind  the 
aircraft. 

In  the.se  figures  the  width  of  the  jet  is  varied  from  0.  25  to  1.  0  times  the 
trailing  edge  span.  The  approximate  lift  distribution  on  the  wing  is  also  plotted 
in  order  to. show  graphically  the  order  of  magnitude  affect  of  the  lift  on  the  flap. 

For  this  purpose  he  thickness  effects  and  the  fuselage  effects,  if  any,  ?re  ignored 
which  results,  using  conical  flow  theory,  in  a  constant  spanwise  averaged  lift 
distribution. 


As  aircraft  Mach  numbers  increase,  the  exhaust  volume  flow  from  conventional 
engines  becomes  larger  due  to  the -decrease  in  air  specific  impulse.  At  Mach  12 
the  mass  flow  from  two  cylindrical  exhausts  would  fill  the  entire  base  of  the 
aircraft.  Engine  cowl  drag  would  become  a  predominant  factor  in  this  case  and 
would  obviously  force  a  compromise  between  cowl  drag  and  expansion  ratio.  The 
lift  distribution  for  aircraft  having  two  cylindrical  exhausts  is  shown  in  Figure  2.4.  16. 
Note  that  at  Mach  12  the  width  of  the  jet  exhaust,  is  equal  to  the  trailing  edge  span 
and  hence  the  je‘_  flap  is  shifting  "lift"  far  behind  the  trailing  edge  of  the  wing. 
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2.4.  3  FINITE  ASPECT  RATIO  EFFECTS 


The  preceding  calculations  have  Leen  conducted  with  the  assumption  that 
the  downwash  behind  the  wing  can  be  neglected  as  in  two-dimensional  theory.  If 
the  downwash  is  taken  into  consideration,  however,  the  exhaust  jets  are  more 
slowly  bent  into  the  direction  of  the  undisturbed  horizontal  flow.  In  frictionless 
flow,  and  with  exhaust  jet  spans  smaller  than  the  trailing  edge  spar,  of  the  wing, 
the  jet  flow  is  completely  imbedded  in  the  downwash  flow  and  will  asymptotically 
reach  not  the  direction  of  the.  undisturbed  flow  but  the  direction  of  the  downwash 
flow  behind  the  wing.  At  very  large  distances,  the  downwash  and  the  jet  will 
eventually  impinge-  upon  the  ground  and  will  be  turned  into  the  horizontal  direction. 

iv 

The  turning  of  the  downwash  flow  at  the  impingement  on  the  ground  will 
pioduce  a  pressure  phenomenon  which  can  be  integrated  as  a  secondary  sonic 
boom  signature  which  moves  along  the  ground  with  the  Macli  number  of  the 
aircraft,  but  with  much  larger  rise  times  than  are  experienced  in  the  primary 
sonic  boom.  (See  Figure  2.4.  17). 

The  slower  turning  of  the  jet  flow  in  the  downwash  has  the  result  that  the 
jet  lift  is  shifted  much  further  behind  the  solid  surfaces  of  the  wing  than  is 
shown  in  figures  2.  4.  1 3  to  2.  4.  16.  Consequently,  the  sonic  boom  alleviation 
effects  are  more  pronounced  when  the  finite  aspect  ratio  effects  are  taken  into 
consideration.  It  is  recommended  to  examine  numerically  the  effects  of  down- 
wash  _>n  the  lift  shift  behind  the  wing  and  the  secondary  sonic  boom  signature  in  the 
impingement  region  of  the  downwash  on  the  ground. 


FIGURE  2.4.1.  2  Dimensional  flat  plate  at  small  angle  of  attack 
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FIGURE  2.4.2.  2  Dimensional  flat  plate  in  ground  effect 
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FIGURE  2.4.3.  2  Dimensional  flat  plate  at  small  angle  of  attack 


FIGURE  2.4.4.  Lift  on  2  dimensional  wings  of  various  chord  lengths 
and  same  lift 


FIGURE  2.4.5. 
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IGURE  2.4.10.  I ift/orag  ratio  and  specific  air  impulse 


FIGURE  2.4.13.  Jet  flap  pressure  distribution  without  downwash 
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FIGURE  2.4.14.  Jet  flap  pressure  distribution  without  downwash 
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expansion  at  20°  jet  deflection) 


FIGURE  2.4.16.  2  Cylindrical  jei  exit  configuration 


CONCLUS  ONS  AND  RECOMMENDATIONS 


3.  1  Considerable  progress  was  made  in  this  contract  in  developing  a  new 

theory  which  permits  the  use  of  normal  size  models  in  wind  tunnel  tests 
to  determine  the  sonic  boom  on  the  ground  by  measuring  the  pressure 
signature  in  the  wind  tunnel  test  section.  The  theory  has  been  developed 
for  axisymmetric  flow'  conditions-  and  a  start  has  been  made  on  the 
theory  for  non-symmetric  lifting  configurations.  It  is  recommended 
that  the  non-symmetric  lifting  theory  be  completed  and  that  procedures 
for  checking  the  theory  in  a  wind  tunnel  be  developed  and  suitable 
verification  tests  conducted. 

3.2  Some  calculations  were  carried  out  by  means  of  an  improved  method 
for  determining  aircraft  contours  with  desirable  finite  pressure  rise 
times;  included  also  was  the  determination  of  phantom  body  shapes 
which  have  predetermined  finite  rise  times  which  can  be  simulated 
by  heat  addition  to  the  flow  upstream  of  the  aircraft,  for  instance  by 
lasers. 

Concerning  the  application  of  laser  beams  to  produce  a  desired  heat 
input  distribution  through  the  mechanism  of  focusing  laser  beams  and 
ionization  of  the  flow,  a  few  experimental  test  points,  from  previous 
tests  with  different  objectives,  are  available.  It  is  recommended  that 
brief  laboratory  tests  be  designed  and  conducted  in  order  to  determine 
to  what  extend  the  control  of  the  heat  input  at  desired  specific  locations 
can  be  accomplished  in  laboratory  and  in  flight  operations. 


3.3  Many  previous  investigations  have  shown  that  deviations  from  the 
eonvemional  circular  engine  exhaust  can  cause  dramatic  improve¬ 
ments  of  the  maximum  lift  coefficient  and  large  effects  on  the  noi^e 
level  o:  engine  exhausts.  The  research  under  this  contract  has  shown 
that  a  -slot -type  engine  exhaust  has  the  unique  capability  of  shifting 
lift  from  the  solid  surfaces  of  the  wing  to  the  region  behind  the  wing, 
as  far  as  the  effects  on  the  external  flow  and  on  the  sonic  boom  are 
concerned.  These  results,  which  for  the  first  time  have  clarified 
an  important  characteristic  of  slot  nozzles,  have  beer,  developed  for 
two-dimensional  flow  conditions.  It  was  shown  that  not  only  a  shifting 
of  the  lift  on  the  main  wing  and  an  associated  reduction  of  the  sonic 
boom  was  achieved,  but  also  favorable  effects  on  the  performance  of 
the  aircraft  in  certain  flight  regimes  and  configurations  were  obtained. 
In  view  of  the  current  interest  in  the  application  of  non-circular  engine 
exhausts  for  maximum  lift  and  noise  attenuation  purposes,  it  is  recom¬ 
mended  that  the  two-dimensional  theory  for  the  slotted  nozzle  exhaust 
effects  be  extended  to  the  three-dimensional  theory,  to  be  used  as  the 
basis  for  subsequent  verification  testing  in  wind  tunnels  and  eventually 


in  flight. 
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